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INCOMPRESSIBLE FLUTTER CHARACTERISTICS OF REPRESENTATIVE AIRCRAFT WINGS 1 

By G. H. Wilts 


SUMMARY 

The present report gives the results oj a detailed study of the 
flutter characteristics of four representative aircraft wings. This 
study was made using the electric-analog computer at the Cali- 
fornia Institute of Technology. During the course of this 
investigation eight important parameters of each wing were 
varied and, in addition, the effects of mass, inertia, pitching 
spring, and location of a concentrated mass were investigated 
for all four wings and at several sweepback angles. 

The introduction of this report discusses in general terms the 
flutter characteristics of airplanes. The second section con- 
tains a discussion of the electric-analog principles that made a 
study of this magnitude feasible. The third section contains a 
discussion of the aerodynamic and structural approximations 
made for simplifying the flutter analysis of a wing. The fourth 
section gives information relating to the errors introduced by 
the finite-difference approximations to continuous aeroelastic 
systems. In addition, data are given pertaining to the flutter 
characteristics of a swept-aoing wind-tunnel model, and the 
results of computations based on two assumptions regarding 
aerodynamic forces on a swept wing are also given. The fifth 
section lists the physical characteristics of the four representa- 
tive aircraft wings and the sixth section contains the computed 
flutter characteristics of the four wings. 

INTRODUCTION 

Flutter i8 a phenomenon which is observed in the transient 
or unforced response of an aerodynamic system. Mathe- 
fnaticaUy speaking, it is observed in the solution of the 
homogeneous differential equation describing the behavior 
of an airplane in flight through still nonturbulent air. An 
airplane wing which is considered to be a continuous beam- 
like or platoliko structure has an infinite number of degrees 
of freedom, and the characteristic equation -which describes 
the transient response has an infinite number of roots. 
Experience has shown that only the roots of lower magnitude 
(frequency) exhibit the problem of instability or flutter. It 
is this fact which makes it possible to predict flutter using 
an analog computer, which represents only the lower fre- 
quency modes of the structure, or using a few normal modes 
in either digital or analog computation. 

The exponents in the transient response of a linear system 
are the roots of the characteristic equation. Since the char r 
acteristic equation involves real parameters, the roots are 
real or occur as complex conjugate pairs. The latter roots 


are the ones of interest here. The real part of a conjugate 
pair is the .reciprocal of the time constant in the transient 
response and the (positive) imaginary part is the frequency 
of oscillation. This is illustrated in figure 1. Mathematical 
description of the transient term is 

y=A x e (a+i ' "> *■ +A 2 e' (<r ~ { u) ' 
or in terms of real functions 

y=Ae°‘ cos (a d+f) 

If the real part of the pair of roots <r is negative the “tran- 
sient” dies out and the root is said to be stable. If the real 
part is positive the transient grows exponentially until 
limited by nonlinearities or destruction, and the root is said 
to flutter. The terminology is not strictly correct, but it is 
common practice to refer to the exponents of the transient 
response as flutter roots since they are numerically equal to 
the roots of the characteristic equation. Throughout this 
report such terminology will be used. 

Damping of flutter roots may bo measured by two dimen- 
sionless numbers f and g, which differ from each other by a 
factor of 2. The former is generally used by control-system 
engineers; the latter, by flutter analysts. They can be de- 
fined by the equation for the particular term in the transient 
response given earlier 

y=Ae° t cos (c ctAd>)=Ae~ ta ^ cos [7(1— f 2 ) w »^+^] 

=Ae 2 cos {V[l— (gV4)]<*»t+<l>} 



Figure 1. — Transient response corresponding to roots of characteristic 

equation. 


■ Supersedes NAOA TN 3780, “Incompressible Flutter Characteristics of Representative Alrcralt WIngB” by O. H. Wilts, 1937. 
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Flutter computations are usually centered around regions 
where the value of g lies in the range. — 0 . 2 <j 7 < 0 . 2 . In 
such cases the factor N /l— J* differs from unity by less than 
0.5 percent. For this reason it is customary to omit this 
factor in the trigonometric term giving the following 
approximation: 

gy 

y=Ae 2 cos (co„t-j-<£) 

This practice will be followed in this report. For damping 
which is small, an approximate rule of thumb is that the 
damping factor g is nearly equal to the per unit decrement 
per cycle divided by t. If percent decrement per cycle 5 is 
used, there results the convenient approximation 

~ 5 

9 IOOtt 

The flutter roots of an airplane are complex functions of 
all geometrical, structural, and inertial properties of the air- 
frame as well as of the airspeed and air density. "With all 
other properties held constant, the lowest airspeed at which 
the flutter root exhibits neutral stability is called the flutter 
speed. If g is plotted as a function of velocity, the abscissa 
(speed) at which the curve first crosses the axis < 7=0 is the 
flutter speed. In this study such curves were used to deter- 
mine the flutter speed, but such curves are used in this 
report only to illustrate the behavior of some unusual flutter 
roots. A tabulation of flutter speeds does not always give 
a good picture of the flutter characteristics. An example is 
shown in figure 2, where the damping of two roots is shown. 
One root becomes unstable at a speed of about 300 miles 
per hour and the other, at a speed of about 600 miles per 
hour. If a parameter variation increases the damping g of 
both roots by 0.03, one flutter speed is raised to 350 miles 
per hour, a 17-percent increase; the other is raised to 603 
miles per hour, a 0.5-percent increase. A further increase 
in g of 0.02 will raise the second flutter speed 0.4 percent to 
605 miles per hour, while the first-root will now exhibit no 
flutter. It should be emphasized that even though a design 
speed of, say, 500 miles per horn: has been surpassed, the 
system may still be regarded as unsatisfactory. A system 
so dose to flutter at a speed of 360 miles per hour might 
actually flutter because of weight (fud) variations or minor 
differences in stiffness resulting from variations within the 
manufacturing tolerances. From the standpoint of this re- 
port all three of the sets of roots discussed above will be 
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Figobb 2. — Variation of g with velocity for typical flutter roots. 


regarded as having essentially the same flutter characteris- 
tics even though they exhibit radically different theoretical 
flutter speeds. Emphasis is given to this point because 
remarks to be made later in this report may be misunder- 
stood without a clear conception of this viewpoint. 

This investigation was conducted at the California Insti- 
tute of Technology under the sponsorship and witli the 
financial assistance of the National Advisory Committee for 
Aeronautics. 


SYMBOLS 


A, Ai, Ai 

Cl. 

d 

E 

El 


constants _ 
half chord 
half chord at root 
half chord at tip 

symbolic representation of circulatory compo- 
nent of lift force due to angle of attack 
lift coefficient 
wing station from root, in. 

Young’s modulus of elasticity 
equivalent beam flexural rigidity, (lb) (sq in.) 
experimental; used as subscript 
flutter frequency, cps 

normal mode frequency of cantilevered engine 
and nacelle, cps 

flutter frequency for continuous structure 
shear modulus 

equivalent beam torsional rigidity, (lb) (sqin.) 

damping factor of a damped sinusoid, e 2 cos oil 

vertical deflection, positive down, in. 

moment of inertia per unit length, lb sec 3 

torsional stiffness 

increase in stiffness, percent 

radius of gyration, in. 

semispan of wing 


// 

U 

u 
0 
QJ 

9 
h 
I 
J 

j 
k 
l 

M, Mi, M it M z , M it M h twisting moment about elastic 

axis per unit length of wing, 
positive nose up, lb 

to mass per unit length, lb sec 3 /sq in. 

m c mass of concentrated mass 

riif fuselage mass 

m w ' total wing mass, lb sec 3 /in. 

mwf total wing mass outside of fuselage 

to lumped mass 

P, Pi, P„ Pi, Pi lift force per unit longth of wing, positive 
nose down, lb /in. 

f Laplace transformation variable 

q n dynamic pressure based on normal component 

of velocity, (1/2) pv 2 , lb/sq in. 
t time, sec 

v airstream velocity, in./sec 

v b flutter velocity of airplane with bare wing 

v f airstream velocity at which flutter occurs, in./sec 

v n component of airstream velocity perpendicular 

to elastic axis, v cos A, in. /sec 
v 0 reference velocity, in./sec 

flutter velocity for continuous wing 
w distance measured along wing 

x 0 distance from midchord aft to elastic axis, in. 
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x t distance from quarter chord aft to elastic axis, 

x„+(b/ 2), in. 

x 3 distance from three-fourths chord forward to 

elastic axis, b —x u in. 

x 3 distance from elastic axis aft to center of 

mass, in. 

y general variable 

A y cell size for finite-difference structure 

a absolute pitch angle about elastic axis, positive 

nose up, radians 

5 percent decrement per cycle 

f per unit critical damping 

6 slope of elastic axis or roll about horizontal 

axis normal to elastic axis, positive tip down, 
radians 

A sweepback angle of elastic axis, deg 

P air density, lb-sec a /in. 4 

r wing twisting gradient, ba/bw 

<r real part of pair of roots 

a angular frequency, radians/sec 

un undamped natural frequency, radians/sec 

(’), (") derivatives with respect to time 

ELECTRIC-ANALOG METHODS OF FLUTTER ANALYSIS 

The use of electrical analogs for the solution of aeroelastic 
problems has been discussed in detail in reference 1. The 
purpose of the present section is to summarize the principles 
briefly. For purposes of flutter analysis the structural 
system is assumed to be linear, and a linear electrical network 
is constructed whose electrical behavior approximates the 
dynamic behavior of the linearized structure. For this 
purpose, capacitors are ordinarily used to represent con- 
centrated or lumped inertia properties, inductors are used 
to represent lumped flexibility properties, and transformers 
are used to represent the geometrical properties of the 
structure (refs. 1 and 2). In such electrical analogs voltages 
throughout the network represent velocities in the structure 
and currents represent forces. Electronic equipment is used 
to produce currents which depend on voltages in the elec- 
trical system in the same manner in which aerodynamic 
forces depend upon the velocities of the airfoil. 

The composite electrical structure can be regarded as an 
electrical model of the aircraft in the same manner that a 
wind-tunnel model would be regarded as a structural model. 
The advantage of this approach lies in the relative ease with 
which one can alter the properties of the model, thus per- 
forming flutter computations with great rapidity. It should 
bo emphasized that the normal modes of the structure are 
not used as tools or elements in the analysis. The analysis 
consists, in fact, in observing the behavior of an electrical 
model of an aircraft in flight. 

That behavior which is most readily observed is the tran- 
sient response to a sudden disturbance. This method is 
therefore similar to the testing technique which is ordinarily 
used for wind-tunnel models. An advantage of the electrical 
method is that tuned pulses may be used so that separation 
of two or more nearly unstable or slightly unstable modes of 
oscillation is more readily accomplished. Basic recorded 
data consist of the logarithmic decrement of the response 
and the frequency of oscillation. Flutter speed and fre- 


quency for any configuration are ordinarily found by com- 
puting the damping g and frequency / for specific values of 
velocity and interpolating to find the frequency and speed 
at which g is zero. 

APPROXIMATIONS FOR SIMPLIFYING FLUTTER ANALYSIS 
STRUCTURAL REPRESENTATION 

For dynamic analysis of airplane wings of large aspect 
ratio, it is customary to treat the wing as a beamlike struc- 
ture in both vertical bending and torsion. It is usually 
assumed for simplicity that an elastic axis exists. For an 
unswept wing this is a straight line which undergoes no ver- 
tical displacement when the wing is subjected to a pure 
torque parallel to this axis and along which no twisting gra- 
dient exists when vertical loads are applied anywhere along 
this line. For an unswept wing of conventional construction 
this simplification is usually quite accurate. For a swept 
wing an elastic axis may be defined as a straight line which 
assumes a constant slope over its entire length when a 
twisting moment is applied parallel to this line and which 
has no twisting gradient when vertical loads are applied 
anywhere along this line. For aspect ratios greater than 5 
or 6 and for conventional wing construction, a line can be 
found on the structure which satisfies this definition reason- 
ably well except near the root. It is not uncommon to find 
an equivalent elastic axis at about the 35- or 40-percent chord, 
a line located aft of the leading edge a distance equal to 35 
or 40 percent of the local chord. 

The assumption of an elastic axis involves the tacit assump- 
tion that chordwise bending of the wing is negligible. It 
follows, then, that the motion of the wing at any spanwise 
coordinate can be described by two coordinates, the vertical 
displacement of some point on the chord, and the angle of 
twist of the chord. If wing motion is described in terms of 
vertical motion of the elastic axis and twisting motion about 
this line, then these motions are not coupled through the 
action of elastic forces in the wing except in the root region 
for a swept wing. 

The root region of a swept wing is necessarily a relatively 
complicated structure. However, for aeroelastic problems 
an equivalent simple structure can be found which is com- 
pletely satisfactory for wings of large aspect ratio. This can 
be demonstrated by the following reasoning. The outer 
sections of a wing exhibit definite beamlike properties, but 
in the region of the root considerable warping of the wing 
surface must take place. The aerodynamic forces near the 
root of the wing are therefore not adequately described by 
strip theory (see discussion of strip theory under “Aero- 
dynamic Forces”). In addition, the inertia effects of this 
section are not readily computed. However, the effects of 
the aerodynamic forces on the root section are insignificant 
for ordinary flutter computations. This has been demon- 
strated many times with the analog computer by removing 
the aerodynamic forces on the inboard cell of the finite- 
difference structure. The inertia forces are also insignificant 
compared with the elastic forces transmitted by the root 
section, and it is therefore possible to replace this section for 
purposes of analysis by a set of “influence coefficients” 
relating transmitted forces to displacement of an outer sec- 
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tion of the wing relative to the fuselage. It has been found 
that in some cases these influence coefficients resemble co- 
efficients for a simple beam extending straight into the 
fuselage and attaching there in some simple way. The wing 
structural axis then consists of a short section which may be 
perpendicular to the fuselage center line and which is simply 
attached to a swept back elastic axis which extends to the 
■wing tip. 

Methods for determining the equivalent structure are 
outside the scope of this report. Since this structure varies 
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(a) Fighter A. (b) Fighter B. (c) Bomber A. (d) Bomber B. 
Figure 3. — Plan forms and cell division of basic airplane. 


greatly with the particular wing construction used, it was 
necessary to choose a simple though typical root structure for 
this study. The structure chosen is illustrated in figure 3 
where the elastic axes are shown by broken lines. The break 
in the elastic axis is assumed to be at the edge of the fuselage, 
and the axis inside the fuselage is assumed to bo straight and 
perpendicular to the airplane center line. The wing is 
assumed to be pinned at the side of the fuselage. Conse- 
quently, all twisting moment is removed at this point and 
it is not necessary to make any assumptions regarding twist- 
ing rigidity inside the fuselage. Bending rigidity inside the 
fuselage is, however, important for symmetric motion. 
During the past 6 years extensive flutter computations have 
been made with the electric-analog computer for commercial 
and military aircraft as well as for wind-tunnel models 
includi n g those described in references 3 and 4. In all coses 
investigated, it has been found that relatively large varia- 
, tions in root conditions have a negligible effect on the flutter 
characteristics (in the sense described in the Introduction). 
Observed changes in damping were usually in the range 
O<|A0[<B.O5, which has very small effect on flutter speed 
unless the curve of g against velocity is very flat, near zero 
values of g. Needless to say, both symmetric and anti- 
symmetric motion of the airplane must be permitted since 
the flutter characteristics for the two types of motion may 
be quite different. 

Fuselage stiffness and inertia properties usually have 
values such that an assumption of a rigid fuselage for sym- 
metric motion alters the flutter characteristics little. For 
fighter planes the error introduced is negligible. For large 
bombers the change in flutter speed may be appreciable, but 
it does not alter the trends to be observed upon variation of 
wing properties. It has therefore been assumed in this 
study that the airplane fuselage is rigid. Tail-surface 
flexibility does not significantly affect wing flutter problems. 
A rigid tail surface with sufficient area to provide satisfactory 
static stability has therefore been assumed. 

AERODYNAMIC FORCES 

For all the flutter computations given in this report, the 
aerodynamic forces have been simplified by two important 
assumptions: 

(1) The air flow is incompressible. 

(2) If the airfoil is divided into strips perpendicular to 
the elastic axis, then the forces on each strip can be com- 
puted as a function of the normal component of the air- 
stream velocity and the motion of that strip can be computed 
independently of the motion of adjacent strips. 

The first assumption is not required by analog methods in 
general, but its use greatly increases the rapidity with which 
data can be obtained. Since the purpose of the study is not 
to obtain specific accurate flutter speeds but to study trends 
in flutter characteristics, this assumption does not seem 
unreasonable. With regard to the use of strip theory two 
assumptions are often found in the literature. In using the 
“airstream method” the wind is divided into strips parallel 
to the airstream, and the forces and moments on each strip 
are computed as though the wing were not swept and the 
air flow about the section were a two-dimensional incom- 
pressible flow. The aerodynamic coefficients may be taken 
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to be the same as those for an unswept wing or may be modi- 
fied by a factor cos A. In applying the “normal-component 
method,” the wing is divided into strips perpendicular to 
the elastic axis. The aerodynamic forces and moments are 
computed as though the effective air velocity were the 
normal component v cos A, and the forces depend only on 
the motion of the individual strip and not upon the motion 
of adjacent strips (except that some small terms may be 
included which are proportional to the twisting gradient and 
therefore dependent upon the motion of the nearest strips) . 
A critical discussion of the two alternatives is given in 
reference 6. This reference recommends use of the normal- 
component method. 

Before adopting the second assumption, an effort was 
made to find some correlation with experimental results. 
Reference 3 contains experimental flutter speeds for a wind- 
tunnel model wing with sweepback angle equal to about 35°. 
This angle is sufficient to give an appreciable difference in 
results obtained with the various assumptions mentioned 
above. The section entitled “Finite-Difference Errors” in 
the present report contains the results of computations which 
show that the normal-component method gives results which 
are as satisfactory as those given by any other method used. 

Equations for determining the aerodynamic forces by this 
method are given in reference 5. In the equations given 
there, several terms are found whose theoretical justification 
is not well established. These terms (grouped in special 
brackets on p. 16 of ref. 5) were found to have negligible 
effect on sample flutter computations. It seems reasonable, 
therefore, to omit these terms from computations involved 
in the present trend study. With these omissions and with 
obvious changes to conform to the symbols and notation 
used in the present report, the equations are: 

P=Pi+P>+P* 

iWa=ilp ~\~M 2 -(-ills 

Pi=-2?r(<Zn) (2&)c(0 [|+(«+« to A)+g (a+6 tan A)] 

P- 2 =-2Hqn) (£) (A) 

Pz=—irpb t (h—X 0 a) 

P 4 =- 2v( 2n ) (2b) [<?(~) (~ 0+x 2 t)+£ 0“^)] A 

M t =— a-iPi 
M a =— xJPi 

M 3 2 ff ( 2 0 (0 (a) 

M*= — irpb* a— aA-J 

=2.( 5 o m [ 2 (a) 1 i+ V )+ 

l T+ S <_ (¥ +I *')3 t “ A 

The terms are grouped in the order shown for convenience in 
establishing analog circuits. The last term in P : is not 
found in the corresponding equation of reference 5. This 


term is removed (mathematically) by insertion of an equal 
but opposite term in P 4 and a similar term in Ms. 

It is added to Pj here because the circuits which generate the 
term a+0 tan A also provide the term (x 2 /v n ) (a+6 tan A), 
the last part of which is not found in reference 5. As is 
pointed out below, this term has a negligible effect so that 
its inclusion is of no importance, but it is indicated in the 
expression for P, for the sake of completeness. It should be 
emphasized that the dynamic pressure g n is based on v„, where 
v n is the velocity component normal to the elastic axis. The 
coordinates a and 9 are both measured in elastic-axis coordi- 
nates. The symbolism C(bp/v n ) is used to represent the 
Theodorsen or Wagner function. A short discussion of the 
interpretation of this symbolic representation can be found 
in reference 6. 

All terms found above can be represented by simple 
analog circuits with the exception of P 4 and M t . Examina- 
tion of equations 6 and 7 of reference 5 shows that each term 
in P 4 and Ms is similar to (if not equal to) a term found in 
the special brackets. Since the latter terms have been 
omitted, there seems to be no logical reason for retaining 
P 4 and M s . Inasmuch as their inclusion greatly complicates 
the analog circuits, these terms were also omitted. 

In addition to the finite-difference approximations and 
those contained in the assumptions of incompressible flow 
and strip theory, three other aerodynamic approximations 
should be mentioned. The first of these is the failure to 
modify aerodynamic forces at the wing tip. The delay in the 
growth of lift forces, as described by the Wagner or Theodor- 
sen functions for two-dimensional flow; cannot apply near 
the tip. Indeed, both the delay in lift and the magnitude 
of the lift must go to zero at the tip. The extent of the error 
introduced depends upon the importance of tip forces in 
flutter computations. Insofar as then- location is concerned, 
these forces are quite important, but, because of wing taper, 
the magnitude of the total force per unit length diminishes 
near the tip. Since wings of considerable taper are involved 
in this investigation, it is to be expected that the error will 
be relatively small. The second approximation is failure to 
compute aerodynamic forces properly at the root of a swept 
wing. As stated earlier, the error introduced by this approxi- 
mation is negligible since the aerodynamic force for a large 
section of the wing root can be omitted entirely without an 
appreciable change in flutter speed. The third approxima- 
tion is introduced by the necessity of computing the Wagner 
function (or the Theodorsen function) electrically. This 
function is computed using networks shown in reference 1 
with an error no greater than 2 percent over the frequency 
range or time interval of interest. 

FINITE-DIFFERENCE ERRORS 

It is the purpose of this section to summarize work which 
was carried out at the Analysis Laboratory of the California 
Institute of Technology to determine finite-difference errors 
in flutter computations for several specific structures. 

FINITE- DIFFERENCE STRUCTURES 

No practical methods have been devised for representing 
general continuous structures with continuous electrical 
systems. The electric-analog computer utilizes lumped elec- 
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trical elements which can, in principle, be used only to con- 
struct analogs for lumped mechanical systems. However, 
as pointed out in references 1 and 2, it is possible to repre- 
sent the dynamic characteristics of beamlike structures by 
a lumped structure based upon finite-difference approxima- 
tions to partial differential equations. It is convenient to 
call this lumped system a finite-difference structure, whether 
it is a mechanical model or an electrical analog. These 
references outline the process by which inertia and stiffness 
properties and aerodynamic forces are averaged or replaced 
by single concentrated inertias, springs, or forces in the 
finite-difference structure. 

It should be remarked, at this point, that the assumption 
of a finite-difference structure insures a finite number of 
flutter roots or exponential functions in the transient re- 
sponse, whereas the continuous structure has, in principle, 
an infinite number. Since the higher frequency roots have 
high damping, only the lower frequency roots are of interest. 
Two or three of these may, however, show essentially zero 
damping simultaneously at a given velocity, and it is some- 
times necessary to determine the characteristics of several 
flutter roots. There is obviously a lower limit to the number 
of cells that must be used to obtain satisfactory accuracy, 
since each cell adds roughly two roots to the system. 

There is little information in the literature which pertains 
to the accuracy with which such structures represent the 
continuous system. Reference 7 gives data for static- 
deflection and normal-mode characteristics of certain finite- 
difference structures but no information about accuracy of 
flutter computations. 

Using equations for aerodynamic forces based on two- 
dimensional strip theory and linear incompressible fluid flow, 
several “exact solutions” have been obtained for flutter 
problems. Some of these are found in references 8 and 9. 
These solutions are exact in the sense that no further physical 
or mathematical simplifications are involved and the only 
errors are introduced by round-off errors in evaluating trans- 
cendental functions and infinite series. Solution of these 
same problems by use of finite-difference approximations to 
partial differential equations provides the most practical way 
of estimating finite-difference errors for other configurations 
for which exact solutions are not obtainable. It is true that, 
in all cases mentioned above, the airfoil has been assumed to 
have uniform spanwise properties and that in most practical 
cases the airfoil has a significant taper. On the other hand, 
reference 7 contains a study of the finite-difference errors in 
the deflection characteristics and normal-mode properties of 
both uniform and tapered beams. This study showed no 
unusual differences in these properties, and so it is assumed 
that the results obtained for flutter of uniform airfoils are 
typical of results that would be obtained for flutter of tapered 
airfoils. 

UNIFORM AIRFOIL WITH PINNED ENDS 

A uniform beam with pinned ends will support ODly sinu- 
soidal modes in both bending and torsion. Flutter modes 
are also of sinusoidal shape and it is therefore possible to 
reduce the flutter problem to an eigenvalue problem which 
can be solved with a high degree of numerical accuracy. 
The finite-difference analogs for a pinned-pinned beam like- 


TABLE I.— PHYSICAL CHARACTERISTICS OF A PINNED- 
PINNED AIRFOIL 


Length, in 288 

Half chord, b, in 48 

Mass per unit length, m, lb secVsq in 0. 01036 

Inertia per unit length, I, lb sec 2 6. 170 

Bending rigidity, El, lb in. 1 1. 412X10* 

Torsional rigidity, OJ, lb in. 1 0. 87X10 9 

*i, in 0. 0 

Xj, in —9. 0 

Air density, p, lb sectyin. 4 0. 0846X10 -8 


wise will support only sinusoidal modes. It is possible 
therefore to get exact solutions for the finite-difference 
approximations to the continuous airfoil. 

The airfoil chosen for this analysis is described in table I. 
For the continuous wing, the flutter speed and frequency 
were found to be C/=692 miles per hour and jf= 12.72 cycles 
per second, respectively. Analysis of the finite-difference 
structure was carried out using eight-, four-, and two-cell 
divisions between the pinned ends. Results are given in 
table II and figure 4. For this particular case it is necessary 
to use more than four cells if flutter speed is to be obtained 
with less thaD 2-percent error. By use of symmetry condi- 
tions at the center of the beam it is necessary to use only 
half this number of cells with an electric-analog computer. 
Thus, use of two analog cells gives a theoretical error of 
about 2.2 percent, and use of four analog cells gives an error 
of ODly 0.6 percent. 




Figure 4. — Finite-difference flutter-speed errors of pinned-pinned 

beam. 
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TABLE II.— COMPUTED FLUTTER SPEED AND FREQUENCY 
OF FINITE-DIFFERENCE PINNED-PINNED AIRFOIL 


Number 
of cells 

Vf, 

mph 

Atv/tico 

//, 

cps 

Af/lf CO 

03 

692 

0 

12. 7 

0 

8 

688 

.006 

12. 6 

. 008 

4 

677 

.022 

12.3 

.032 

2 

644 

. 069 

11. 0 

. 134 


UNIFORM CANTILEVER WING WITH CONCENTRATED MASS 

Analytical determination of the flutter speed of a cantilever 
wing is much more difficult than that for a beam with pinned 
ends. However, other investigators have obtained accurate 
numerical solutions for a few configurations. The most 
important of these is described in reference 9. This case is 
of importance for two reasons : It involves several spanwise 
positions of a large eccentric concentrated mass which greatly 
affects the flutter speed; and for some positions at least two 
completely different flutter roots can be found. 

Table III presents the physical characteristics of the airfoil 
analyzed in reference 9. In this reference, the flutter speed 

TABLE III— PHYSICAL CHARACTERISTICS OF UNIFORM 
CANTILEVER WING 


Half chord, 6, in 4 

Span, l, in - 48 

Mass per unit length, m, lb secVsq in L 877 X10 -4 

Pitching inertia per unit length, lb sec 5 8. OOXIO -4 

Flexural rigidity, El, (lb) (sq in.) 1. 407 X10 5 

Torsional rigidity, OJ, (lb) (sq in.) 0. 692 X 10 5 

Elastic-axis position, in —0. 604 

Center-of-mass position, Xj, in 0. 158 

Mnss of concentrated mass, m c , lb sec’/in a 23 X 10~® 

Pitch inertia of concentrated mass, lb secVin 0. 1636 

Center of mass of concentrated mass, (x 3 ) „ in —3. 274 

Air density, p, lb sec 1 /in. 4 1. 165 X10 -7 


arid flutter frequency were computed for seven mass locations, 
data for which are reproduced in table IV. Since the location 
of a concentrated mass may be important in flutter analysis, 
and since all points on a finite-difference beam are not 
equally suitable as an attachment point for a concentrated 
mass, it was believed that a comparison of the above data 
with finite-difference solutions was quite important. Un- 
fortunately, similar accurate solutions for a finite-difference 
structure are not readily obtained, so it was necessary to use 
the electric-analog computer to obtain these solutions. The 
resulting comparison therefore contains both finite-difference 
and analog-computer errors. Previous work has indicated 
that the latter are probably not greater than 1 percent if the 
Theodorsen function is represented accurately. 

In this analysis two slightly different beam analogs were 
used. In both, the beamlike properties were represented by 
a system of levers (transformers), but in one group the lumped 
forces were applied at the junctions of the levers and in the 
second group the forces were applied at the midpoints of the 
levers. The analog of the second group was once thought 
to give a better approximation since it resembles the Russell 
beam analog discussed in reference 7. Recent investigation 
has shown that this belief is without foundation, and the 
second analog is now preferred only as a matter of conven- 
ience for sweptback wings since it provides the wing slope 
directly at the force stations where it is needed for computa- 
tion of aerodynamic forces. In both cases the cantilever 
condition at the root was provided by a half cell at the root, 
and the forces nearest the tip were applied one-half cell from 
the tip. Thus the first group involved an integral number of 
cells, and the second group involved a half integral (integer 
plus one-half) number of cells. Five cases were investi- 
gated; 2, 2 i, 4, 5i, and 6 cells. Since it was shown in the 
previous section that less than 4 cells was of no interest for 
present purposes, only the results of 4, 5§, and 6 cells are 
presented in this report. 


TABLE IV.— THEORETICAL AND EXPERIMENTAL FLUTTER CHARAC- 
TERISTICS OF CANTILEVER WING WITH CONCENTRATED MASS 


[Data taken from references 9 and 11] 


Mass lo- 

Calculated 

Experimental 

cation ; 







per unit 

Vf, 

Vf, 

fr, 

Vf, 

Vf, 

//> 

span 

mph 

per unit 

cps 

mph 

per unit, 

cps 



C) 



(*) 


0 

227 

0. 400 

25. 27 

228 

0. 401 

1 

. 167 




221 

. 388 

19. 1 

. 229 

226 

. 397 

19. 23 

221 

. 388 

17. 4 

. 292 




233 

. 410 

16. 3 

. 333 




266 

. 451 

15. 5 

. 354 

277 

. 488 

2R 04 

260 

. 458 

b 16. 3-26. 8 


o 359 

. 631 

30. 68 



( d ) 

. 938 

273 

.481 

25. 67 

261 

.469 

( d ) 

. 959 

251 

. 442 

24 87 

261 

. 442 

2L 8 

. 979 




231 

.407 

2L 6 

1. 000 

206 

. 360 

23. 60 

218 

. 384 

2L 4 


» 1.0 per unit velocity Is 1CH im/sec or 668 mph. 

a This experimental record seems to show nearly simultaneous divergence and flutter at two frequencies, 
c Calculated divergence speed Is about 279 mph. However, a flutter speed can still be calculated mathematically, 
a Divergence was observed experimentally. 


020507 — 00 90 
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In view of the simplicity of the flutter curves shown in 
reference 9, it was expected that data would be taken at 
only a few spanwise mass locations. However, it was soon 
found that the flutter characteristics were much more com- 
plicated than anticipated, and data were taken at 24 mass 
locations in the 6-cell case. The flutter characteristics of 
the 'wing with variable location of the concentrated mass 
are sketched in figure 5 (a). As the concentrated mass is 
moved outward from the root, the flutter speed drops slightly. 
At a distance about 16 percent of the total span from the 
root a m i nim um is reached, and beyond the 25-percent posi- 
tion the flutter speed rises very rapidly. At the 30-percent 
position the flutter speed for this root has become equal to 
the flutter speed of a completely different root. The flutter 
speed for this second root drops with increasing spanwise 
position of the mass making it impossible to determine with 
the analog computer the speed for the original root beyond 
the 30-percent position. The flutter speed for the second 
root reaches a minimum with the mass at the 45-percent 
position, then rises to a very high value as the mass is moved 
toward the 75-percent position. A flutter root which is 
probably the second is observed for mass positions near the 
tip, the lowest flutter speed occurring with mass at the tip. 


It was also observed that divergence of the wing occurred 
whenever the flutter speed exceeded about 5,000 inches per 
second. Because of divergence, it was not possible to meas- 
ure with accuracy flutter speeds which exceeded divergence 
speed by more than about 50 percent. As a result, flutter 
speeds with mass near the 75-percent span could not be 
measured. 

The flutter characteristics for the 4-, 5£-, and 6-cell struc- 
tures are shown in table V and figure 5 (b). Data for the 
seven positions analyzed in reference 9 are also plotted in 
the figure. Inspection of these curves shows that many more 
accurate numerical solutions are required to determine the 
finite-difference errors for all mass positions. In spite of 
the inadequate numerical data, an attempt was made to draw 
a smooth curve through the known points taken from refer- 
ence 9. In doing this the 5£- and 6-cell analog data were 
used as a guide in determining the shape of the curve. This 
curve, shown in figure 5(a), has already been discussed. It 
is realized that a significant error of as much as 2 or 3 per- 
cent may exist in this curve for some mass positions, but there 
was no other method for obtaining estimated errors for the 
finite-difference structures. With the understanding that 
the comparison data may be in error in some regions, figure 6 



(a) Theoretical and wind-tunnel flutter characteristics. 

(b) Flutter characteristics of finite-difference analog. 

Figttbb 5. — Flutter characteristics of uniform cantilever wing with concentrated mass. 
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TABLE y,— FLUTTER CHARACTERISTICS OF FINITE- 
DIFFERENCE ANALOG OF UNIFORM CANTILEVER 
WING 


Mass position 

Vf, mph 

vt, per unit 
(•>) 

//. cps 

4 cell 

0 

224 

0. 394 

23. 9 

. 126 

212 

. 373 



. 160 

213 

. 374 



. 226 

225 

. 396 

18. 9 

. 260 

241 

. 424 

la i 

. 276 

275 

. 484 



.325 

269 

. 473 

27.7 

. 360 

268 

. 454 

27. 2 

. 376 

246 

. 431 

26. 3 

.600 

273 

. 480 

27. 6 

. 600 

297 

. 522 

2R 2 

. 626 

303 

. 633 

2a 4 

. 660 

376 

. 661 

3L 4 

.876 

327 

. 676 

26. 7 

. 960 

237 

. 417 

24, 0 

1. 000 

197 

. 346 

22. 4 

5# cell 

0 

223 

0. 392 

23. 8 

. 182 

218 

.383 

20. 3 

. 236 

227 

. 399 

ia 6 

. 264 

236 

. 415 

17.8 

. 273 

249 

. 439 

17. 3 

. 291 

296 

. 520 

2a5 

.309 

289 

. 508 

2 a 0 

. 328 

284 

. 600 

2a 3 

. 364 

274 

.482 

27. 9 

.464 

244 

. 430 

27. 2 

. 646 

284 

. 600 

2a 3 

. 600 

313 

. 560 

29. 7 

. 636 

333 

. 686 

30. 7 

. 856 

377 

. 664 

29. 4 

.909 

305 

. 637 

27. 0 

.946 

258 

. 464 

26. 6 

L 000 

212 

. 373 

2 a 9 

6 cell 

0 

228 

0. 401 

2a 7 

.083 

223 

.393 

23. 3 

. 167 

219 

.386 

20. 8 

. 217 

226 

. 398 

ia 7 

. 250 

239 

. 420 

ia o 

. 267 

266 

. 448 

17. 1 

. 284 

286 

. 604 

16. 7 

. 317 

296 

. 521 

2 a 2 

. 367 

270 

.476 

27. 6 

.417 

252 

.443 

26. 8 

. 460 

268 

. 464 

27. 2 

.600 

266 

. 468 

27. 6 

.550 

277 

. 487 

27. 8 

.583 

282 

. 496 

2 a 2 

. 600 

306 

. 639 

2a9 

. 616 

334 

. 588 

30. 1 

. 667 

»477 

«.840 



.868 

349 

. 614 

2 a 4 

. 883 

319 

. 562 

27. 4 

.917 

276 

.486 

2a 3 

. 960 

242 

. 426 

25. 1 

. 965 

226 

. 398 

24, 6 

.984 

213 

. 376 

24. 2 

L 000 

201 

. 353 

23. 6 


• Moss position Is per unit span measured from root 
b 1.0 per unit velocity- 6GS mph. 



Figure 6. — Flutter-speed errors of finite-difference analog of uniform 
cantilever wing with concentrated mass. 

was prepared showing the percentage error in flutter speed 
for the various analogs as functions of the mass location. 
For 5|- and 6-cell structures the average errors are about 2 
percent. It can he readily seen that, although a 4-cell 
analog gives very satisfactory results for the hare wing (mass 
position 0), it is necessary to use more than 4 cells if errors 
less than 5 percent are required at other mass locations. A 
further discussion of this investigation will he found in 
reference 10. 

As a result of this analysis, it was decided that all flutter 
computations made in this trend study would he made using 
6| cells to represent one-half of the airplane wing. 

EXPERIMENTAL CORRELATION 

"Wind-tunnel tests have been made of many model struc- 
tures. It is difficult, however, to find unclassified data in 
which the structure is completely and accurately described. 
In the course of this investigation, two cases were found in 
which a correlation between experimental and computed 
characteristics could be attempted. The first of these is the 
uniform unswept cantilever wing discussed in the preceding 
section. The flutter speed and frequency observed in a 
wind tunnel are reported in reference 9 and a companion 
report, reference 11. These data are summarized in table 
IV, which also contains the computed values of reference 9. 
A better understanding of the correlation is obtained if the 
experimental data are plotted with the assumed analytic 
solution. Figure 5 (a) shows such a comparison. The corre- 
lation for this case seems unusually good. 




1394 


REPORT 1390 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


PLOTTER SPEED OF A SWEPT- WING MODEL 

The second case for which correlation with experiment is 
possible is found in reference 3. This reference gives results 
of wind-tunnel tests to determine the flutter speed of a model 
wing with sweepback angle equal to 34.5°. This wing had 
two concentrated masses attached at approximately the 30- 
and 80-percent span positions. In an effort to compare the 
airstream and normal-component aerodynamics for flutter 
computations, an electrical analog was constructed for this 
wing. For any sweepback angle it is to be expected that the 
two methods will give flutter speeds differing by a factor of 
approximately (cos A) w , unless the aerodynamic coefficients 
are modified by the factor cos A in the airstream method, in 
which case the two methods should give similar results. The 
principal difficulty encountered was determination of the 
properties of the concentrated masses on the wing, since 
reference 3 does not give complete information about these 
masses and their geometrical location. The best data that 
could be deduced from this report are given in table VI. 
Since the masses are alined with the airstream but are repre- 
sented in elastic-axis coordinates, a product of inertia be- 
tween roll and pitch exists. Since no such information was 
available, the product of inertia was omitted from computa- 
tions, and the rolling inertia about a chord line was assumed 
to be one-half as large as the pitching inertia about the 
elastic axis. It is believed that these approximations and 
simplifications will affect the results by less than 1 percent. 

A comparison of observed and computed characteristics is 
given in table VH. The first three normal-mode frequencies 
show satisfactory agreement, with differences of 1, 5, and 3 
percent, respectively. The flutter speed computed with 
either representation of aerodynamic forces is lower than the 
wind-tunnel value. In the case of the airstream method, 
the discrepancy is 19 percent, or, if the aerodynamic coeffici- 
ents are modified, 11 percent. Using the normal-component 
method, the discrepancy is 12 percent. Flutter frequency is 
in error about 20 percent in all cases. Although the observed 
differences are relatively large in all cases, it is concluded 
that the normal-component method recommended in refer- 
ence 5 is satisfactory for this model. 

CHARACTERISTICS OF FOUR REPRESENTATIVE AIRCRAFT 

Plan forms and stiffness and inertia data were chosen after 
surveying the various fighter, bomber, and transport planes 
developed in recent years. Four representative airplanes 
were chosen, two fighters and two large bombers. Smaller 
attack bombers and transports were not included because 
of lack of time. The airplanes chosen are not similar in all 
respects to any particular set of four airplanes, but they do 
have stiffness and inertia properties which resemble four 
specific aircraft. Plan form, sweepback angle, elastic-axis 
location, and concentrated-mass locations were, however, 
chosen more arbitrarily so that this report could remain 
unclassified. The four basic plan forms are shown in figure 3. 
The basic fighter A has a bare unswept wing with span of 
about 500 inches, taper ratio of 2.0, and aspect ratio 6. The 
basic fighter B has a wing sweepback angle of 30°, a span 
of about 400 inches, and a taper ratio of 2.0. The two basic 
wings have the same length measured along the elastic axis 


TABLE VI.— PHYSICAL CHARACTERISTICS OF SWEPTBACK 
WING WITH CONCENTRATED MASS 
[Data taken from reference 3; more detailed information will be found in this reference. Mass 
of wing is for portion outboard of root restraint. Data for concentrated masses are not 
given explicitly in reference 3 and must be regarded as only approximate.] 


Wing characteristics: 

•Span, in 48. 3 

b Root half chord, b„ in 5. 2 

b Tip half chord, b x , in 2. 30 

Wing total mass, m., lb secVin 0. 00784 

Tunnel fluid density, p, lb secVin.* 3.40X10 -7 

Sweepback angle, A, deg 34. 5 

Concentrated-mass characteristics: inboard Outboard 

Mass, m c> lb secVin 0. 00806 0. 00462 

Pitch inertia about elastic axis, lb sec s /in._ 0. 0712 0. 0192 

Per unit spanwise position (from root) 0. 30 0. 78 

Center-of-mass position, (rj) n in — 1. 74 0. 60 

* Measured along elastic axis. 
b Measured perpendicular to elastic ails. 


TABLE VH.— EXPERIMENTAL AND COMPUTED FLUTTER 
CHARACTERISTICS OF SWEPTBACK WING "WITH CON- 
CENTRATED MASS 


Type of frequency 

Normal mode frequencies, 
ops, at — 

Mode I 

Mode 2 

Mode 3 

Experimental model frequency 

6. 97 

30. 9 

37. 9 

Measured analog frequency 

6. 91 

32. 6 

39. 1 


Flutter characteristics 

Type of result 

v/, mph 

A Vf/Vf, 

//, cps 

Wind-tunnel results, r/,_ 

193 

0 

20. 1 

Analog results, normal-component 
method 

170 

-. 12 

24. 2 

Analog results, airstream method *„ 

157 

-. 19 

24 0 

Analog results, airstream method b ._ 

173 

-. 11 

24 0 


* Lift coefficient Coo— 3r. 
b Lift coefficient Clo—Ix cos A. 


and the same chords measured perpendicular to the elastic 
axis. 

The basic bomber A has an unswept wing with span of 
about 1,700 inches, taper ratio of 2.5, and aspect ratio 12. 
It has a concentrated mass representing an engine nacelle 
at the 0.46-span position with center of mass about one-half 
chord forward -of the elastic axis. The basic bomber B has 
a wing sweepback angle of 30°, a span of about 1,500 inches, 
and a taper ratio of 2.4. It also has a concentrated mass 
representing an engine nacelle at the same relative position 
as for bomber A. The two basic wings have the same length 
measured along the elastic axis and approximately equal 
chords when measured parallel to the airstream. 

Mass per unit length, pitch inertia per unit length, bend- 
ing rigidity, and torsional rigidity were drawn as smooth 
curves approximating the characteristics of some typical 
modern aircraft. As described in reference 2, these data must 
be collected or lumped over distances corresponding to the 
cell length of the analog finite-difference structure. The 
assumed curves and the lumped values are shown in figure 7. 
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The l um ped values are also listed in tables AM L 1 to XI, which 
give all pertinent characteristics of the basic airplanes. 

Eight important parameters of the basic airplane wings 
were varied in an effort to find similar features in the flutter 
characteristics of the various wings. The quantities varied 
and the extent of their variation are summarized as follows: 


Minimum Maximum 

Qnantlty varied value value 

Wing mass density, per unit basic 0. 5 2. 0 

Wing pitch inertia, per unit basic 0. 5 2. 0 

Bending rigidity, per unit basic 0. 67 1. 5 

Torsional rigidity, per unit basic 0. 67 1. 5 

Center-of-mass location, percent chord 25 60 

Elastic-axis location, percent chord 30 50 

Chord, per unit basic 0. 67 L 5 

Sweepback angle, deg 0 45 


With the exception of sweepback angle, these quantities 
were varied one at a time from their basic value. However, 
for all four basic airplanes, some or all of the parameters 
were varied for two or three values of sweepback angle. It 
is realized that the above variations do not constitute a 
comprehensive survey. However, to a considerable extent 
the changes in flutter speed due to several variations are 
additive if the variations are small and are made simultane- 
ously. Another limitation is that the flutter characteristics 
are affected by the spanwise variation in the first seven 

TABLE Vni.— CHARACTERISTICS OF BASIC FIGHTER A 
(n) Physical characteristics 


Sweepback angle, A, deg 0 

Semispan of wing,* l, in 238 

Cell size for finite-difference structure, Ay, in 34 

^Root chord, b 6„ in 106 

Tip chord, b b x , in 63 

Taper ratio 2- 00 

Aspect ratio 6 

Wing elastic axis, percent chord 40 

Wing center of mass, percent chord 40 

Total wing mass external of fuselage, m wp lb sec’/in 10. 7 

Fuselage mass, tny, lb seoVin 21 

Fuselage radius of gyration, pitch, 0 in 100 

Fuselage radius of gyration, roll, in 25 

Fuselage center of mass aft of elastic axis, in 0 

Tail center of pressure aft of elastic axis, in 230 

Tail area, sq in 3, 000 

Air density, p, lb secVin. 4 1. 146 X10 -7 


(b) Inertia and stiffness values lumped for finite-difference structure 


Station 

1 

2 

3 

4 

6 

6 

Per unit span * 

Half chord, b 6, in 

Lumped mass, m, lb 

secVin 

Lumped pitch inertia 

0. 214 
47. 3 

1. 73 
695 

0. 367 
43. 5 

1.34 

568 

1 

0. 643 
36. 0 

0. 66 
253 

0. 786 
32. 2 

0. 40 
126 

0. 928 
28. 4 

0. 23 
53 

nQ">f(,dy/EI) 

<>10 '°f(dylGJ) 

19. 8 
56.3 

61.3 

145 

112 

219 

205 

365 

415 

667 


850 


» Measured along elastic axis. 

•> Measured perpendicular to elastic axis, 
e About elastlo axis. 

d Stiffness values are lumped between mass stations. 


TABLE IX.— CHARACTERISTICS OF BASIC FIGHTER B 
(a) Physical characteristics 


Sweepback angle, A, deg 30 

Semispan of wing,* l, in 238 

Cell size for finite-difference structure, Ay, in 34 

Root chord, b br, in 100 

Tip chord, b b x , in 63 

Taper ratio 2. 00 

Wing elastic axis, percent chord 40 

Wing center of mass, percent chord 40 

Total wing mass external of fuselage, lb sec J /in 14. 00 

Fuselage mass, m/, lb sec^/in 21 

Fuselage radius of gyration, pitch, 0 in 100 

Fuselage radius of gyration, roll, in 26 

Fuselage center of mass aft of elastic axis, in 0 

Tail center of pressure aft of elastic axis, in 230 

Tail area, sq in 3, 000 

Air density, p, lb secVin. 4 1. 140X10 - ’ 


(b) Inertia and stiffness values lumped for finite-difference structure 


Station 

1 

2 

3 

4 

6 

0 

Per unit span *_ _ 

Half chord, b b, in 

Lumped mass, m, lb 

secVin. . — 

Lumped pitch inertia 

0. 214 
47. 3 

1. 62 
628 

0. 357 
43.6 

1. 43 
577 

0. 600 
30. 8 

1. 32 
603 

0. 643 
36. 0 

1. 16 
411 

0. 786 
32. 2 

0. 94 
310 

0. 928 
28. 4 

0. 06 
204 

n0'°f(dylEI) 

d 10 '“/(dy/GM)- 

33.0 
71. 4 

100 

193 

156 

283 

221 

661 

321 

400 

447 

430 

224 


* Measured along elastlo axis. 

Measured perpendicular to elastlo axis. 

• About clastic axis. 

d Stiffness values aro lamped between mass stations. 

quantities listed. The four basic cases, two fighter and two 
bomber airplanes, can be regarded as four oases in which 
spanwise variations have been made. However, since five or 
six of the quantities are varied in going from one case to 
another, it is not possible to determine the effect of span- 
wise variation of only one of the parameters. Other quanti- 
ties which were thought to have second-order effeots were 
not considered. Am ong these are altitude (represented by 
ratio of air density to wing mass), fuselage mass and pitching 
inertia, and tail configuration. This does not imply that 
flutter velocity is independent of altitude, but with very 
minor variations the flutter velocity varies inversely as the 
square root of the air density. Sea-level air density was 
used throughout this study. 

It is improbable that bombers of the plan form and size 
studied will be flown without engines on the wing. Conse- 
quently, the basic cases of interest are those in which a 
concentrated mass is located there. On the other hand, it is 
of some interest to compare the characteristics of the baro 
wing as well as those of a wing with concentrated mass. 
Both bombers A and B were studied with bare wing as woll 
as with concentrated mass in the basic position on the wing. 

For purpose of reference it is necessary to assign a number 
to designate each particular case. The group discussed above 
comprises 175 cases. The assignment of case numbers is 
shown in table XII. This table shows most rapidly the 
various cases that were studied. 
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TABLE X.— CHARACTERISTICS OF BASIC BOMBER A 
(a) Physical characteristics 


Sweepback angle, A, deg 0 

Somispan of wing,* l, in 845 

Cell size for finite-difference structure, Ay, in 130 

Root chord, b b n in 200 

Tip chord, b b z , in 80 

Taper ratio 2. 50 

Aspect ratio 12 

Wing elastic axis, percent chord 40 

Wing center of mass, percent chord 40 

Total wing mass external of fuselage, triwp lb secVin 39. 7 

Fuselage mass, m/, lb secVin 120 

Fuselage radius of gyration — pitch, 0 in 240 

Fuselage radius of gyration — roll, in 50 

Fuselage center of mass aft of elastio axis, in 0 

Tail center of pressure aft of elastic axis, in 700 

Tail area, sq in 20, 000 

Air density, p, lb secVin. 4 1.146 X 10 -7 


(b) Inertia and stiffness values lumped for finite-difference structure 


Station 

1 

2 

3 

4 

5 

6 

Per unit span,* . . 

0. 154 

0. 308 

0. 461 

0. 616 

0. 769 

0. 923 

Half chord, b 6, in 

90. 8 

81. 5 

72. 3 

63. 1 

53. 9 

44 6 

Lumped mass, m, lb 








secVin 

6. 28 

6. 46 

4 16 

2.47 

0. 98 

0. 60 

Lumpod pitch intertia. 

13, 300 

10, 600 
1 

7, 150 
1 

3, 640 
1 — 

1, 170 

360 

1 


M0 

*W°f{dylGJ)_.J 


3. 3 
10.0 


& 4| 
23. 1 


14.3 

30.9 


30. 5 
47. 6 


101 

143 


650 

618 


618 


• Measured along elastic ails. 
i> Measured perpendicular to elastio ails. 

0 About elastic axl3. 

a Stillness values arc lumped between moss stations. 


Concentrated masses on fighter wings usually consist of 
fuel tanks, bombs, or similar stores. It is impossible, there- 
fore, to select a single value for mass and inertia which can 
bo regarded as typical. For certain positions many values 
for mass and inertia were chosen, although in most of the 
cases studied the number of values was restricted by the 
time available for computations. For reference purposes, 
the basic mass for fighter planes was arbitrarily chosen to 
bo one-quarter of the mass of the entire wing (half of the mass 
of one side), the pitching radius of gyration was set equal to 
30 inches, and the roll radius of gyration was assumed to be 
15 inches or less. Specific data for the two fighters are listed 
in table XIII. 

Concentrated masses for bomber airplanes are usually en- 
gine nacelles, with a mass which can be predicted within a 
factor of 2. Nevertheless, it is of some interest to study the 
effect of various mass values in these cases also. Basic mass 
value for both bombers was assumed to be 15 pounds-seconds 
squared per inch, which corresponds to a weight of nearly 
6,000 pounds. Pitching radius of gyration was assumed to 
be 35 inches. Basic mass position was assumed to he at the 
0.46-span position and 60 inches in front of the elastic axis . 
Those data are also tabulated in table XIII. 

The concentrated-mass characteristics varied in this study 
are: 


(1) Mass 

(2) Pitching inertia about center of mass 

(3) Spanwise location 

(4) Chordwise location 

(5) Pitching flexibility 

The assignment of case numbers is more difficult for this 
phase of the study. Although specific spanwise positions 
were chosen, it was not possible to choose chordwise positions 
beforehand. The chordwise positions were chosen as the 
data were obtained. In some cases more than 20 positions 
were used for a given spanwise location. Consequently, one 
case number was assigned to all chordwise variations at a 
given spanwise location. A summary of all variations with 
the corresponding case numbers is given in table X I V. 

Pitching flexibility of the concentrated mass was varied 
in six cases involving both bombers. In all cases, the chord- 
wise location of the center of mass was basic (table XUI). 
In three cases the mass was in basic spanwise position and 

TABLE XI.— CHARACTERISTICS OF BASIC BOMBER B 
(a) Physical characteristics 


Sweepback angle, A, deg 30 

Semispan of wing,* l, in 845 

Cell size for finite-difference structure, Ay, in 130 

Root chord, b b r , in 170 

Tip chord, b b„ in 70 

Taper ratio 2. 43 

Wing elastic axis, percent chord 40 

Wing center of mass, percent chord 40 

Total wing mass external of fuselage, m^ p lb secVin 48. 7 

Fuselage mass, mj, lb seo^/in 120 

Fuselage radius of gyration, pitch,” in 240 

Fuselage radius of gyration, roll, in 60 

Fuselage center of mass aft of elastic axis, in 150 

Tail center of pressure aft of elastio axis, in 700 

Tail area, sq in 20, 000 

Air density, p, lb secVin. 4 1.146 X 10 -7 


(b) Inertia and stiffness values lumped for finite difference structure 


Station 

1 

2 

3 

4 

5 

6 

Per unit span*., 

Half ohord, b b, in . 

Lumped mass, m, lb 

secVin. 

Lumped pitoh inertia 

0. 154 
77. 3 

6. 28 
3, 940 

0. 308 
69. 6 

5. 23 
3, 080 

0. 461 
61. 9 

4 26 
2,200 

0. 615 
54 2 

a 36 

1, 560 

0. 769 
46.5 

2. 74 
910 


MO' 0 f(dy/ El) 

M0 ">f(.dy/GJ) 

7. 12 
6. 31 

18. 8 
15. 5 

31. 0 
22. 6 

64 9 
42.2 

91. 6 
77. 8 

|J|m 

■ 


• Measured along elastio ails. 

• Measured perpendicular to elastic nils. 

• About elastic ails. 

a Stillness values are lamped between mass stations. 

in three cases the mass was at the tip. Case numbers are 
given in table XIV. 

TRENDS IN FLUTTER CHARACTERISTICS 

REFERENCE QUANTITIES AND GRAPHICAL PRESENTATION 

Results of the study of trends in flutter characteristics 
which are listed in table XV are given in miles per hour 
















TABLE X3L— ASSIGNMENT OP CASE NUMBERS 




Case number of — 

Variable factor 

Valuo or type 
1 

Fighter 
A ; 

A-0° 

Fighter 

A = 45° 

Fighter 

A=0 ,J 

Fighter 

A=i5° 

Bomber 

A; 

A=0° 

*M 

Bomber 

A; 

A=30°; 

»M 

Bomber 

A; 

A— 0^; 

»B 

Bomber 

B -' „ 
A=30°; 

•M 

Bomber 

A=6«; 

Bomber 

B; 

A-=45°; 

*M 

Bomber 

A-^Sb°; 

*B 

gig 

Bomber 

B; 

A-46°; 

•B 

Basic 


1 


23 

42 

58 

72 

84 

08 

mm 



152 

163 

Mass factor 

2 

2 


24 

43 

69 

74 

85 

69 

mm 

IKP5R 


164 

106 


1/2 

3 


26 

44 

00 

75 

80 

100 

■Si 

127 


166 


Fitch, inertia factor 

2 

4 


26 

45 

61 

76 

87 

10 i 

mm 



i i f ft 
100 

167 


1/2 

6 


27 

• 40 

02 

77 

88 

102 

■ifi 

129 

141 

167 

108 

Bending rigidity factor 

3/2 

6 


28 

47 

03 

78 

89 

103 

119 

130 

142 

158 

109 


2/3 

7 


29 

48 

04 

79 

90 

104 

120 

131 

143 

169 

■779: 

Torsional rigidity faotor 

3/2 

8 


30 

49 

05 

80 

91 

106 

121 

132 

144 

160 

171 


2/3 

9 


31 

32 

60 

00 

81 

92 

100 

122 

133 

145 

101 

172 

Looai stiffness variation 


10 


67 

93 

b Center of mass, percent- 

25 

32.5 

40.. 


20 

19 


61 

62 

42 




— 

— 

— 

— 

— 

— 


1 

23 

58 


84 

98 

113 

114 

187 

162 

163 


46. 

11 

33 

63 

08 

94 

107 

123 

134 

146 

162 

173 


50 

12 

21 

34 

64 

69 

82 

96 

108 

124 

135 

147 

103 

174 


00. 

13 

22 

35 

65 

70 

83 

96 

109 

125 

136 

148 

104 

175 

b Elnfltio ads, percent 

1 - 



30 

50 












1 


23 

42 

58 


84 

98 



137 




50 

14 


37 

67 

71 


97 

110 



149 



Chord factor. . 


15 


38 


111 



160 




2/3 

10 


39 



Rife* 


112 

Riitali 

| ■ jj | R 

151 


^■R 

Sweepbaolt angle, deg 

o’..„ 

1 


23 

23 

68 

68 


113 

113 

113 

152 

162 

152 


15 ... 

17 


40 

40 







IS 


41 

41 

72 

72 


98 

98 

98 

137 

137 

137 


45. 

19 


42 

42 

73 

73 


114 

114 

114 

163 

163 

163 

Nacelle pitch flexibility 

Maes at station 3.. 



176 



178 

177 






Mass at tip 

— 

— 

— 

— 

179 

— 

— 

181 

180 

— 

— 

— 

— 


* M, concentrated nissi ociricg; B, btra wing. 

b Center of maa tod tlastie arte mtararod hi percent chord from leading edge. 
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TABLE Xin.— CONCENTRATED-MASS CHARACTERISTICS 
AND LOCATIONS 

(a) Characteristics 



Fighter 

Bomber 

Basic mass, lb sec 3 /in 

2. 68 

15. 0 

Pitch radius of gyration,* in 

30 

35 


15 

17. 5 


6 

7 

Roll radius of gyration, in 

15 



0 


Basic spanwiso position, 11 in _ 


390 

Basio chordwiso position, 0 in 

— 

60 


* About clastic axis; alrstream coordinates. 

■> Outboard from center line measured along elastic axis. 
« Forward of clastic axis, parallel to alrstream. 


(b) Location 


Airplane 

Locations of concentrated mass, per unit span,* at 
station — 

1 

2 

3 

4 

4.6 

5 

6.5 

6 

Tip 

Fighter. 

Bomber. 

0. 214 
. 164 

0.368 
. 308 

0. 500 
. 461 

0. 643 
. 615 

a 692 

0. 786 
. 769 

a 846 

0.929 
. 923 

L 00 
1. 00 


» Distances are measured In per unit span along elastlo axis. 


and in per unit values of a reference speed. The reference 
velocity chosen is 

d o = 10 4 in./sec=568 mph 

Obviously a flutter speed of 1.5 would not represent a real- 
istic value since this would correspond to supersonic speed 
with a Mach number of about 1.1. However, such a number 
still has useful significance for two reasons: (1) A major 
purpose of this study is to establish trends and to determine 
what configurations tend to be more or less stable than 
others, and (2) a change in stiffness is equivalent to a change 
in velocity, so that a structure with one-half the stiffness of 
another, but otherwise unchanged, would exhibit a flutter 
speed 1/V2 times as great as that of the other, a value equal 
to 1 .06 or 600 miles per hour in the case given above. 

All geometrical, structural, and inertia quantities are given 
in per unit values. For example, distances are measured in 
units of the airplane semispan, and masses in terms of a basic 
value. Specific numerical values can be obtained using 
figure 3 which gives linear dimensions of the airplanes, figure 
7 which gives inertia per unit length and rigidity values, 
and tables VJLH to XI which list all other pertinent charac- 
teristics of the four basic airplanes. The density of air at 
sea level was used throughout these computations. The 
value chosen is : , 

p=1.146X10- 7 lb sec 2 in." 4 


In presenting results graphically, flutter speeds have, in 
general, been reduced to dimensionless values by using the 
flutter speed of the basic configuration as the velocity unit. 
For example, when plotting antisymmetric flutter speed hs 
a function of wing mass density for a particular wing such 
as that of fighter B with A=45°, the flutter speeds have 
been divided by the antisymmetric flutter speed of fighter 
B, A=45°, with basic wing mass. Symmetric and antisym- 
metric results are both presented, rather than choosing the 
one which gives lowest flutter speed. Where such results 
are presented in the same figure, symmetric results are gen- 
erally indicated by solid lines, and antisymmetric results, 
by dotted lines. Specific numerical values for the flutter 
speeds and flutter frequencies are found in table XV. 

MASS AND INERTIA VARIATIONS 

In most practical configurations the normal mode of vibra- 
tion with lowest frequency is predominantly a bending mode 
and is usually called the first wing bending mode. In the 
absence of a large concentrated mass on the wing, a pre- 
dominant torsional motion is usually observed in the third 
or fourth mode. Simple flutter can often be predicted with 
engineering accuracy using only these two modes as the 
normal coordinates of the structure. When a large con- 
centrated mass is involved, the situation is much more 
complex. Two or more torsion modes as well as two or 
more bending modes become important in flutter compu- 
tations, and several flutter roots may be observed which 
predominantly involve various ones of these modes. For 
eccentric masses it becomes, in fact, impossible to speak of 
bending and torsion modes since many modes will involve 
both largo bending and torsion displacements. 

' In those cases in which flutter involves a bending mode 
and a higher frequency torsion mode, it cau be said that a 
structural change which separates the frequencies of these 
modes ordinarily raises the flutter speed, and a change 
which makes the frequencies more nearly equal lowers the 
flutter speed. This generalization is not always valid. A 
change in mass density without change in pitching inertia 
has greatest effect on first bending frequency even in cases 
with large sweepback. Consequently, increase in wing 
mass density would be expected to give an increase in flutter 
speed, and a decrease in mass density, a decrease in flutter 
speed. Changes in pitching inertia .vould normally be 
expected to have an opposite effect. Such variations were 
made for three fighter configurations, four bare-wing bomber 
configurations, and five bomber configurations with concen- 
trated mass. The mass density and pitching inertia were 
separately changed by factors of 2.0 and 0.5, making a total 
of 48 configurations in addition to the 12 basic cases. Refer- 
ence case numbers are given in table XU. 

Tabulation of flutter speed and frequency for each case 
will be found in table XV. The results are also shown in 
figure 8. As mentioned earlier, the flutter speeds have been 
reduced to dim ensionless values by using as the velocity unit 
the flutter speed of the basic wing for each basic configuration. 
The trends predicted above are found in most cases. In the 
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TABLE XIV.— CASE NUMBERS FOR CONCENTRATED-MASS VARIATIONS 


\ 

*»/»» 6 a.U 

Station 

Fighter A, A=0°, with radius of 
gyration, in., of — 

Fighter A, A =45°, with radius of 
gyration, in., of — 

Bomber A, 
A=0°, with 
radius of gyra- 
tion, in., of — 

Bomber B, 
A=0°, with 
radius of gyra- 
tion, in., of — 

Bomber B. 
A=0°, with 
radius of gyra- 
tion, in., of — 

6 

15 

30 

6 

15 

30 

35 

35 

36 

1 








248 




1 









276 


2 


182 

183 

218 

219 

220 

249 

264 

277 


3 


184 

185 

221 

222 

223 

250 

266 

278 


4 




224 

225 

226 

251 

266 

279 


4. 5 







252 




5 

186 

187 

188 

227 

228 

229 

253 

267 

280 


5 5 







254 




0 

189 

190 

191 

230 

231 

232 

255 

268 

281 


Tip 

192 

193 

194 

233 

234 

236 

266 

269 

282 

0. 60 

1 









283 

9 







257 

270 

284 


3 







258 

271 

286 


4 







259 

272 

286 


5 







260 

273 

287 


5. 5 







261 




6 



* 




262 

274 

288 


Tip 







263 

275 

280 

2 

' Tip 

105 

196 

197 







1 

102 

193 

194 

233 

234 

235 




fiO 


lOfl 

199 

200 

236 

237 

238 




. 10 


201 

202 

203 






08 


204 

205 

206 







. 037 


207 

208 

209 







. 016 


210 








.008 


211 


















1 

5 

186 

187 

188 

227 

228 

229 




. 50 




239 

240 

241 




. 16 


212 

213 

■mmm 







. 037 


215 

216 

mam 










wMMm. 







i 

Tip 

* 242 

* 243 

* 244 

* 245 

* 246 

* 247 





Case number for fighter B. 
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(a) Fighters A and B. 

(b) Bomber A. 

(o) Bomber B. 

(d) Bomber B; concentrated mass on wing. 
Figure 8. — Effect of wing mass and inertia variations. 


case of b are-wing fighters the effect is very systematic. The 
average of all cases is given as follows: 


Mass value 

Pltoh Inertia value 

Average change 
v/ f percent 

0.6 

L0 

-12 

2. 0 

1.0 

6 

L 0 

. 5 

8 

1.0 

2.0 

-9 


The effect, though uniform, is quite small. 

The results for bombers show much less consistency. 
For cases both with and without concentrated masses, the 
effect of wing-mass density variation is unpredictable. 
Nearly half of the cases show trends which are opposite to 
that predicted above. The addition of a concentrated mass 
at the 0.46-span position reversed the trend in several cases. 
On the other hand, change in wing pitching inertia did show 
a systematic trend for all bomber cases. On the average a 
change in pitching inertia by a factor of 2 changed the 
flutter speed about 7 percent. 

The following conclusions can be drawn: 

(1) A change of wing pitching inertia shorn a systematic 

trend for all wings although ihe effect is small. 

(2) A change of wing mass shorn a definite trend for 

typical fighters although the effect is small. 

(3) Change of wing mass for typical large bombers 

with or without concentrated masses shows no 

systematic trend. 

STIFFNESS VARIATIONS 

It has been pointed out (e. g., ref. 1, p. 783) that when 
incompressible fluid flow is assumed, a change of stiffness is 
equivalent to a change of velocity insofar as transient 
response of an airfoil is concerned. Consequently, it can be 
said that a uniform increase in stiffness will raise the flutter 
speed by the square root of the factor by which stiffness is 
increased. In most airplanes it is found that the increase 
in torsional rigidity is primarily responsible for the increase 
in flutter speed and that, in general, a change in bending 
rigidity over rather wide limits does not change the flutter 
speed significantly. 

As shown in table XTT , 12 configurations were studied to 
support this conclusion. Since both bending rigidity and 
torsional rigidity were separately changed by factors of 0.67 
and 1.60, there are a total of 48 case numbers assigned to 
this group. The results of this study are listed in table XV 
and presented graphically in figure 9. For ease of compari- 
son flutter speeds are converted to dimensionless values, and 
flutter characteristics for changes in bending and torsional 
rigidity are plotted side by side. In general, it was found 
that change in torsional rigidity by a factor of % or % increased 
or decreased the flutter speed by 20 percent and that a 
similar change in bending rigidity had a negligible effect 
upon the flutter speed. Among the 12 configurations studied, 
the following exceptions to this trend were noted: 

(1) Bomber A, A=0°: In the antisymmetric case both 
bending and torsional rigidity had roughly equal effect, 
flutter speed changing ±10 percent for a rigidity change of 
% or %. 
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TABLE XV.— FLUTTER SPEED AND FREQUENCY 


Case 

Airplane 

Symmetric 

Antisymmetric 

v/, mph 

v f , per 
unit 

Sr, cps 

Vf, mph 

v f , per 
unit 

Sr, cps 

I 

Fighter A 

989 

1. 74 

14 2 

1, 091 

L 92 

a? 

2 


1 , 080 

L<10 

10. 4 

1, 102 

1. 94 

a 13 

3 


892 

L 57 

16. 8 

1 , 068 

L 88 

11. 3 

4 


835 

1. 47 

12 4 

1 , Oil 

1. 78 

a i 

5 


1, 102 

L 94 

15. 0 

1, 142 

2 01 

a 7 

6 


972 

L 71 

14 8 

1, 085 

1. 91 

9. 1 

7 


1, 000 

1. 76 

13. 7 

1, 119 

1. 97 

a 3 

8 


1, 222 

2. 15 

16. 8 

1, 369 

2 41 

10. 2 

9 


795 

1. 40 

12 0 

886 

1. 56 

7. 4 

11 


813 

L 43 

14 0 

989 

L 74 

8. 6 






f 636 

L 12 

16. 6 

12 


710 

1. 25 

13. 9 

1 892 

1. 67 

a 6 






i 546 

. 96 

la o 

13 


619 

1. 09 

12 9 

{ 778 

1. 37 

7. 9 

14 


773 

L 36 

14 1 

858 

L 51 

a 2 

15 


795 

1.4 

12 0 

790 

1. 39 

7. 1 

16 


1, 295 

2. 28 

15. 7 




17 


’ 938 

L 65 

14 4 

733 

1. 29 

17. 2 

18 


949 

L 67 

15. 0 

750 

1. 32 

17. 7 

19 


1, 028 

L 81 

15. 7 

881 

1. 56 

17.6 

20 


1, 415 

2.49 

14 6 

1, 267 

2 23 

16. 1 

21 


847 

L 49 

14 9 

784 

1. 38 

17.0 

22 


790 

1. 39 

14 6 

796 

1. 40 

15. 9 

23 

Fighter B 


1. 89 

7.4 

1, 017 

1. 79 

6. 9 

24 


1, 006 

1. 77 

5. 3 

1, 057 

1. 86 

4. 8 

25 


876 

1. 54 

10.8 

938 

1. 65 

9. 2 

26 


989 

1.74 

a 6 

926 

1. 63 

6. 1 

27 


1, 131 

1. 99 

ao 

1, 074 

1. 89 

7.4 

28 



1. 80 

a 2 

1, 017 

1. 79 

7. 2 

29 


1,040 

L 83 

6.8 

I, 017 

1. 79 

a 8 

30 


1, 273 

2. 24 

a4 

1, 244 

2 19 

a 3 

31 


835 

1. 47 

a 7 

830 

1. 46 

a 9 

33 


881 

L 55 

a o 

898 

1. 68 

a 5 

34 


767 

L 35 

a i 

824 

1. 45 

a i 

35 


653 

L 15 

7.3 

710 

L 25 


36 


1, 244 

2. 19 

7. 8 

1, 227 

2 16 

7.3 

37 


852 

1. 50 

WBM 

813 

1. 43 

a 6 

38 


784 

L 38 


727 

L 28 

• a 5 

39 


1, 358 

2 39 

a 5 

1, 392 

2 45 

a 3 

40 


943 

1. 66 

a i 

760 

1. 32 

11. 2 

41 


932 

L 64 

9. 0 

801 

1. 41 

11. 4 

42 


983 

1. 73 

9.2 

1, 023 

1. 80 

11. 5 

43 


1, 051 

1. 85 

a 7 

1, 131 

L 99 

9. 2 

44 


847 

1. 49 

11. 8 

847 

1. 49 

13. 9 

45 


864 

1. 52 

a o 

1, 080 

1. 90 

10. 0 

46 


1, 074 

L 89 

9. 8 

1, 159 

2 04 

12 6 

47 


983 

L 73 

10. 2 

983 

L 73 

12 9 

48 


960 

1. 69 

a 7 

1, 074 

1. 89 

11. 2 

49 


1, 176 

2 07 

1L 7 

1, 313 

2 31 

ia 8 

50 


801 

L 41 

a 3 

801 

1. 41 

10. 6 

51 


2, 171 

2 82 

11/5 - 

--- -2, 057 

3. 62 

13. 8 

52 


1, 188 

2 09 

a 8 

1, 142 

2 01 

12 0 

53 


898 

1. 58 

9. 3 

1, 017 

L 79 

11. 2 

54 


835 

L 47 

a 9 

1, 307 

2 30 

7. 9 

55 


756 

L 33 

a4 

989 

L 74 

a4 

56 


1, 347 

2 37 

10. 2 

1, 260 

2 20 

12 4 

57 


784 

L 38 

a o 

949 

1. 67 

1L 4 

























Case 

Airplane 

Symmetric 

Antisymmetric 

Vf, mph 

vj, per 
unit 

ft, cps 

V/, mph 

v f , per 
unit 

ft, ops 

68 

66 

60 

61 

62 

63 

64 
66 
66 
68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 

82 

83 

84 
86 
86 

87 

88 

89 

90 

91 

92 
94 
96 

96 

97 

Bomber A 

710 

625 

943 

602 

784 

688 

733 

898 

663 

625 

586 

540 

534 

1, 085 

1, 193 

784 
1, 165 

986 

1, 017 
1, 034 
886 
1, 086 
847 

750 

580 

818 
864 
761 
676 
909 
790 
830 ' 
1, 017 
648 
705 
636 
540 
614 

1. 25 
1. 10 
1. 66 
1. 06 
1. 38 
1. 21 
1. 29 
1. 58 
. 99 
1. 10 
1. 03 
. 95 
. 94 

1. 91 

2. 10 

1. 38 

2. 05 

1. 70 

1. 79 
1. 82 
1. 56 
1. 91 
1. 49 

1. 32 

1. 02 

1. 44 
1. 52 
1. 34 
1. 19 
1. 60 
1. 39 
1. 46 
1. 79 
1. 14 
1. 24 
1. 12 
. 95 
1. 08 

9. 4 
7. 5 
12. 2 
& 2 
10. 2 
9. 7 
9. 3 
11. 4 

7. 9 

8. 8 
& 4 

7. 8 
& 6 

17. 9 

17. 7 

8. 6 

20. 9 

14 5 

21. 8 

17. 4 
19. 3 
23. 7 
14 2 

18. 3 

17. 5 

9.4 
7. 8 
11. 4 

a 0 
10. 8 
9. 6 
9. 6 
11. 8 
7. 8 
9. 5 
a 9 
a i 

9. 8 

869 
648 
972 
682 
1, 040 
761 
1,000 
1, 227 
619 
676 
563 
455 
528 
f 699 

1 835 

| 886 

1 i, on 

663 
903 
f 597 

1 705 

778 
744 
710 
869 
608 
\ 642 

1 676 

1 563 

\ 608 

886 
767 
983 
864 
852 
989 
773 
949 
807 
648 
563 
466 
563 

1. 53 
1. 14 
1. 71 
1. 20 
1. 83 
1. 34 
1. 76 
2 16 
1. 09 
L 19 
. 99 
. 80 
. 93 
1. 23 
1. 47 
1. 56 
1. 78 
1. 15 
1. 59 
1. 05 
1. 24 
1. 37 
1. 31 
1. 25 
1. 53 
1. 07 
1. 13 
1. 19 
. 99 
1. 07 
1. 56 
1. 35 
1. 73 
1. 62 
1. 50 
1. 74 
1. 36 
L 67 
1.42 
1. 14 
. 99 
. 82 
. 99 

as 

7. 0 

ia 8 

7. 5 
9. 9 

a 6 
12 6 

7. 0 

ao 

7. 4 

6. 7 

7. 7 
9. 6 
(*) 

10. 6 
ia 9 

7. 5 
20. 1 
a 2 

15. 4 
9. 9 
9. 7 
9. 0 
11.0 
7. 9 
a 7 
(*) 

17. 4 

a o 

9. 9 

a 7 

9. 6 
7. 1 
12 0 
a i 
10.1 
12. 6 
6. 6 
9.8 
10. 5 
10. 5 
10. 7 

98 

Bomber B 

1, 295 

2. 28 

mm 

1, 261 

2 22 

9. 8 

99 


1, 290 

2. 27 

WBm 

1, 250 

2. 20 

7. 0 

100 



2. 25 

mam 

1, 341 

2 36 

14 3 

101 


1, 216 

2. 14 

Tael® 

1, 176 

2 07 

9. 5 

102 


1, 341 

2. 36 

1G3B 

1, 318 

2 32 

10. 3 

103 


1, 307 

2. 30 

‘1 

1, 244 

2 19 

10. 2 

104 



2. 24 

10. 5 

1, 278 

2 25 

10. 6 

106 


1, 563 

2. 75 

12. 9 

1, 568 

2 76 

ia 0 

106 


1, 068 

1. 88 

8. 9 

1, 017 

1. 79 

a 3 

107 


1, 148 

2. 02 

10. 1 

1, 108 

1. 95 

9.4 

108 


1, 040 

1.83 

9. 2 

1, Oil 

1. 78 

a 9 

109 


898 

1. 58 

a 1 

876 

1. 54 

a o 

110 



1. 80 

9. 3 

1, 000 

L 76 

9. 1 

111 


966 

1. 70 

a 8 

909 

4 60 

a2 

112 


1, 699 

2. 99 

12. 5 

1,750 

3. 08 

12 6 

113 


1, 295 

2.28 

10. 2 

1, 352 

2 38 

as 

114 


1, 489 

2. 62 

11. 7 

1, 432 

2 52 

11. 3 

116 


1,284 

2. 26 

7. 1 

1, 261 

2 22 

6. 6 
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TABLE XV.— FLUTTER SPEED AND FREQUENCY— Concluded 


Case 

Airplane 

Symmetric 

Antisymmetric 

Vf, mph 

v f , per 
unit 

ft, cps 

Vf, mph 

Vf, per 
unit 

ft, cps 

116 

Bomber B 

1,284 

2. 26 

la 

7 

eg 

a 

37 

15. 3 

117 


1,227 

2. 16 

9. 

4 

sm 

a 

18 

8. 4 

118 


1, 341 

2. 36 

10. 

8 

1, 477 

a 

60 

| 9. 6 
jll. 4 

119 


1,284 

2. 26 

10. 

2 

1, 284 

a 

26 

9. 3 

120 


1,250 

a 20 

10. 

6 

1, 347 

a 

37 

11. 3 

121 


1, 534 

a 70 

ia 

0 

1, 653 

a 

91 

13.8 

122 



1. 85 

a 

3 

1,051 

i. 

85 

7. 6 

123 


1, 108 

1. 95 

9. 

2 

1, 074 

l. 

89 

& 3 

124 


983 

1. 73 

a 

6 

955 

i. 

68 

7. 8 

126 


813 

1. 43 

7. 

4 

773 

l. 

36 

7. 0 

126 


1, 523 

a 68 

8. 

8 

1,477 

a 

60 

7. 8 

127 


1, 472 

a 59 

ia 

5 

1, 409 

a 

48 

14. 8 

128 


1, 392 

a 45 

10. 

4 

1, 335 

a 

35 

10. 6 

129 


; 1,557 

a 74 

la 

2 

1, 500 

a 

64 

ia 0 

130 


1, 628 

a 69 

n. 

5 

1, 477 

a 

60 

11. 0 

131 


1, 420 

a so 

li. 

3 

1, 386 

a 

44 

10. 7 

132 


1, 738 

3. 06 

la 

6 

1, 699 

a 

99 

13. 1 

133 


1, 260 

a 20 

9. 

4 

1, 205 

a 

12 

9.0 

134 


1, 364 

2. 40 

li. 

0 

1, 284 

a 

26 

11.0 

135 


1, 261 

a 22 

10. 

2 

1, 188 

a 

09 

10. 6 

136 


1, 119 

1. 97 

- a 

7 

1, 080 

i. 

90 

10. 2 

137 


1, 165 

a 05 

mm 


' 1,244 

a 

19 

10. 3 

138 


1, 023 

1. 80 



1, 239 

2. 

18 

7.0 

139 


1, 261 

a 22 

13. 

7 

1, 206 

a 

12 

14. 2 

140 


1, 108 

1. 95 

7. 

0 

1, 170 

a 

06 

9. 8 

141 


1, 182 

a 08 

7. 

3 

1, 284 

a 

26 

10.4 

142 


1,206 

a 12 

a 

1 

1,267 

a 

23 

10. 7 

143 


1, 108 

1. 96 

6. 

1 

1, 239 

a 

18 

9. 6 

144 


1, 358 

a 39 

7. 

5 

1, 617 

a 

67 

11. 8 

146 


983 

1. 73 

a 

6 

1,034 

l. 

82 

8. 7 

146 


1, 074 

1. 89 



1, 131 

l. 

99 

10. 2 

147 


1, 023 

1. 80 



1,034 

l. 

82 

9. 6 • 

148 


898 

1. 58 



898 

l. 

68 

8. 6 

149 


972 

1. 71 



1, 023 

l. 

80 

9. 7 

150 


881 

1. 55 



909 

l. 

60 

8. 6 

151 


1, 483 

a 61 



1, 665 

a 

93 

10. 8 

162 


1, 216 

a 14 

10. 

1 

1, 216 

a 

14 

10. 4 

153 


1, 114 

1. 96 

6. 

6 

1, 421 

a 

50 

10. 4 

164 


1, 193 

a io 

6. 

5 

1, 216 

a 

14 

6. 7 

155 


1, 159 

a 04 

ia 

8 

1, 170 

a 

06 

13. 9 

156 


1, 136 

a oo 

9. 

4 

1, 142 

a 

01 

9. 2 

157 


1, 261 

a 22 

10. 

5 

1, 261 

a 

22 

10.4 

158 


1, 199 

a n 

10. 

0 

1, 193 

a 

10 

9. 9 

159 


1, 199 

a li 

9. 

9 

1, 182 

a 

08 

10. 1 

160 


1, 466 

a 58 

la 

1 

1, 449 

a 

55 

ia4 

161 


977 

L 72 

a 

1 

977 

i. 

72 

8. 1 

162 


1, 057 

1. 86 

9. 

9 

1, 057 

i. 

86 

9. 8 

163 


949 

1. 67 

9. 

4 

960 

l. 

69 

9. 6 

164 


818 

1. 44 

a 

7 

818 

l. 

44 

a 6 

165 


1, 040 

1. 83 

4. 

3 

1, 438 

a 

53 

7. 16 

166 


1,244 

a 19 

10. 

2 

1, 376 

a 

42 

14 0 

167 


1, 080 

1. 90 

6. 

4 

1, 347 

a 

37 

10. 0 

168 


1, 125 

1. 98 

a 

6 

1, 460 

a 

67 

10. 5 

169 


1, 205 

a 12 

7. 

4 

1, 455 

a 

56 

11. 2 

170 


1, 000 

1. 76 

a 

5 

1, 386 

a 

44 

9. 7 

171 


1, 227 

a 16 

6. 

7 

1, 699 

a 

99 

11. 9 

172 


983 

1. 73 

a 

0 

1, 188 

a 

09 

9. 1 

173 


1, 040 

1. 83 

a 

5 

1, 313 

a 

31 

10. 9 

174 


1, 000 

L 76 

a 

X> 

1, 216 

a 

14 

11. 0 

175 


920 

1. 62 

a 

5 

1, 085 

l. 

91 

10. 6 
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(c) Bomber B; bare wing. 

(d) Bomber B; concentrated mass on wing. 

a function of wing rigidity. 


(b) Bomber A; bare wing and concentrated mass on wing. 

Figvee 9. — Flutter speed as 
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(2) Bomber B, A=45°: Symmetric case same as case (1) 
above. 

(3) Bomber A, concentrated mass at 0.46 span, A=0°: 
In the antisy mm etric case torsional rigidity had a 50 per- 
cent greater effect (± 30-percent change in flutter speed) 
and bending rigidity had a negative effect (-T- 10-percent 
change in flutter speed). 

(4) Bomber A, concentrated mass, A=30°: In the sym- 
metric case, the trend was normal only for decrease in tor- 
sional rigidity and increase in bending rigidity. 

These exceptions do not constitute a major deviation, and 
the trend is considered well established. 

LOCAL STIFFNESS VARIATIONS 

It is not to be expected that the effect discussed in the 
preceding section will be observed if torsional rigidity is 
changed at various stations along the wing. In the absence 
of a concentrated tip mass, any effect on flutter speed must 
vanish for stations near the tip, and presumably the largest 
effect will be observed for stations near the fuselage. Be- 
cause of the great ease with which these data could be 
obtained, the effect of local stiffness variation was obtained 
for several configurations. 

The analog computer requires lumping or averaging of 
inertia and stiffness properties. Consequently, it is possible 
to determine readily only the effect of a stiffness variation 
which must be assumed to extend over the entire length of 
a cell in the finite-difference structure. The basic data 
consist therefore of step curves. To obtain an approximate 
value for the per unit change in flutter speed per unit change 
in stiffness per unit length at any point along the wing, it is 
necessary to draw a smooth curve passing through the step 
curve such that the areas under the two curves are approxi- 
mately equal. It is believed more suitable to present the 
step curve and let the reader do any smoothing his applica- 
tion requires. The configurations studied are listed below: 

(1) Fighter A, A=0°, symmetric and antisymmetric, 

case 10 

(2) Fighter B, A=0°, symmetric and antisymmetric, 

case 32 

(3) Bomber A, A=0°, bare wing, symmetric, case 93 

(4) Bomber A, A=0°, concentrated mass at 0.46 span, 

symmetric, case 67 

Results of this study are presented in figure 10. The 
abscissa of a curve is the spanwise station at which the bend- 
ing or torsional rigidity variation is made. The ordinate is 
the per unit change in flutter speed per unit change in stiffness 
per unit length along the wing. If, for example, the stiffness 
is increased j percent over a distance w along a wing of 
semispan l between the stations d—(w/ 2) and d+(w/ 2), then 
the ordinate of the smoothed curve at the abscissa d/l when 
multiplied by j(vj /l) will give the approximate percent change 
in flutter speed. 

Data for the two fighters show great similarity both for 
symmetric and antisymmetric conditions. The greatest 
effect is obtained by changing torsional rigidity near the 
midspan position, slightly outboard for fighter A and slightly 
inboard for fighter B. The effect of a bending-rigidity 
change was found to be small at all stations. In most cases 



(a) Fighter A. (o) Bomber A; bare wing. 

(b) Fighter B. (d) Bomber A; concentrated mass at 0.40 span. 

Figure 10. — Effect of local stiffness variation on flutter speed. 
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a small negative effect was observed, the flutter speed drop- 
ping slightly as the bending rigidity was increased. 

Bomber A without concentrated mass showed a similar 
trend with the following exceptions: 

(1) Maximum improvement was obtained by changing 
torsional rigidity near the root of the wing (0.25-span 
position). 

(2) Increase in bending rigidity was observed to decrease 
the flutter speed by an amount which was 5 to 10 times 
greater than that for fighter A. 

Addition of a concentrated mass at spanwise station 0.46 
has a great effect on this characteristic. The mass chosen is 
typical for an airplane engine and is sufficiently large so that 
the wing is, to a certain extent, pinned at this point for the 
particular flutter root involved. Consequently, stiffness 
changes inboard of the engine have a negligible effect, and 
changes outboard have an effect very similar to that observed 
for a bare wing of reduced length. 

It should be remarked at this point that the result dis- 
cussed above is not to be regarded as a trend for all con- 
figurations. When the flutter is primarily an outer-wing 
bending-torsion flutter, this result is to be expected. Ex- 
perience has show, however, that occasionally an inner- 
panel torsion mode is involved in flutter, and change in 
torsional rigidity outboard of the nacelle has no significant 
effect. It is unfortunate that such a configuration was not 
investigated for this report. 

CENTEU-OF-MASS LOCATION 

The location of the wing center of mass has a great effect 
upon flutter speed of an airplane wing. The general trend 
is that flutter speed increases as the center of mass moves 
forward. It is not generally true that it is at a constant 
chord location at all spanwise stations. However, for pur- 
poses of studying tho trends, it is necessary to assume some 
basic position for the center of mass. Past experience has 
show that a center-of-mass location near the elastic axis 
(usually slightly aft) is both realistic and typical. For this 
reason the basic position of the center of mass was assumed 
to be the elastic axis or 40 percent chord. Variation in 
center-of-mass location was between the 25- and 60-percent- 
chord points. Thirteen configurations were studied, the 
various center-of-mass locations comprising 53 cases listed 
in table XII. 

The results are listed in table XV and are shown graphi- 
cally in dimensionless form in figure 11. The general trend 
is that the flutter speed increases as the center of mass moves 
forward and decreases as the center of mass moves aft, 
except for center-of-mass locations far behind the elastic 
axis. For positions near the elastic axis the flutter speed 
changes about 3 percent for a shift in center of mass equal 
to 1 percent of the chord. For the extreme aft positions 
(60 percent chord) most of the curves become quite flat, 
and in about four cases the flutter speed has started to rise 
slightly as the center of mass is moved farther aft. On the 
other hand, the curves become very steep for center-of-mass 
locations forward of the elastic axis. In most cases the 
increase in flutter speed was so great that data could not 
be obtained for the 25- and 32.5-percent-chord locations 
because the flutter speed greatly exceeded the divergence 



(a) Fighters A and B. (b) Bombers A and B. 

Figvbe 11. — Flutter speed as function of center-of-mass position in 

percent chord. 


speed. The average percentage change in flutter speed for 
a shift in center of mass equal to 1 percent of the chord 
depends upon location of the center of mass as indicated 
below: 


Center-of-mass location, percent chord 40 50 60 

Change in flutter speed, percent 3. 1 1. 7 0. 8 


One unusual case was noted. The results for fighter B, 
A=45°, in figure 11 (a) show an unusual behavior for aft 
center-of-mass location in the antisymmetric case. A study 
of the frequency of oscillation for each position tends to 
support the conclusion that two different flutter roots are 
involved. In any case the results are anomalous and could 
bear further investigation. 
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v 


(a) Fighter A; A=0°; antisymmetric. 

(b) Bomber A ; A = 30° ; concentrated mass on wing; antisymmetric. 

HF, high frequency. 

Figure 12. — Plot of damping factor g against velocity for some unu- 
sual flutter roots. 

Fighter A, A=0°, shows another unusual characteristic in 
the antisymmetric case. One flutter root disappears as the 
center of mass is moved forward of the 46-percent-chord 
location. This result, shown in figure 11 (a), is more easily 
understood by referring to figure 12 where the curves of g 
against v are plotted for this configuration. A similar case 
shown in figure 11 (b) has two readily observable flutter 
roots, one with low and the other with high flutter frequency. 
Data for both cases are given in figure 11 (b). It is true 
that only the one with the lower flutter speed is of practical 
interest, but for purposes of studying trends both are equally 
important. A plot of g against v for this case is also shown 
in figure 12. 

ELASTIC- AXIS LOCATION 

A main component of the aerodynamic pressures on an 
airfoil is equivalent to a force applied at the quarter chord. 
Consequently, the elastic-axis location relative to the quarter 
chord determines the nature of the coupling between aero- 
dynamic forces and the structure. If elastic-axis location 
alone were changed, both center of pressure and center of 
mass would change with respect to the assumed structural 
a xis . In order to separate the effects due to these two 
changes, the center of mass was moved with the elastic axis 
in the configurations discussed here. Aerodynamic coupling 
in which the center of pressure (quarter chord) is forward 
of the elastic axis may have a destabilizing influence while 
a center of pressure aft of the elastic axis generally has a 
stabilizing effect. 

Elastic-axis locations between the 30 and 50 percent chord 
were used in the 16 cases listed in table XU. Results are 



Figure 13. — Flutter speed as a function of elastio-axis position In 

percent chord. 


given in table XV and figure 13. In all cases the expected 
trend was observed. For an elastic axis near the 40 percent 
chord, the flutter speed changed 3.2 percent on the average 
for a shift in elastic axis equal to 1 percent of the chord. 
This effect is not linear over a wide range; however, the 
flutter speed increases more rapidly as the quarter chord is 
approached and decreases more slowly as the elastic axis is 
moved aft. For an elastic axis at the 50 percent chord, the 
corresponding change in flutter speed was only 1.8 percent. 

CHORD VARIATIONS 

A change in chord of a wing is usually accompanied by 
significant changes in mass, inertia, and stiffness as well as 
changes in other characteristics. In an effort to assess the 
effect of aerodynamic pressures alone," variations were mado 
in which mass, inertia, and stiffness were held constant while 
the chord length was changed. Location of the elastic axis 
was maintained at a constant per unit chord station so that 
the distance between quarter chord and elastic axis changed 
in proportion to the change in the chord length. Since the 
magnitude of the aerodynamic force increases with chord 
length and since the predominantly destabilizing lag of the 
Theodorsen function increases with chord length, it is to bo 
expected that the flutter speed will decrease as the chord 
length is increased. 
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I I I I I I 

5 1.0 15 .5 ID 15 


Per unit basic chord length 

Figure 14. — Flutter speed as function of chord length. 

Four configurations were studied in which the chord 
length was changed by factors of 0.67 and 1.50. The eight 
eases and the configurations are listed in table XII. Flutter 
characteristics are given in table XY and figure 14. The 
results are remarkably uniform. On the average, a 7-percent 
change in flutter speed results from a 10-percent change in 
chord, smaller chords giving a higher flutter speed. 



SWEEPBACK 

The effect of sweepback upon flutter speed depends upon 
many factors. In conventional wing design, the root 
structure varies greatly with sweepback angle, and the 
equivalent elastic axis may show considerable variation in 
position. For wings of large sweepback angle and low 
aspect ratio, the concept of an elastic axis may not be useful 
in describing structural properties. From another point of 
view the problem is even more perplexing since there is not 
general agreement about the nature of the aerodynamic 
forces on a swept wing. In the section of this report entitled 
"Finite-Difference Errors” the results of three methods of 
computation were compared with results of wind-tunnel 
tests of a model wing which was sweptback 34.5°. Two 
methods were found to give similar results, which were 
significantly better than those of the third. Although the 
agreement was not entirely satisfactory, it was decided to 
use the aerodynamic forces recommended in reference 5. 
For the present investigation, the following assumptions 
were therefore made: v 

(1) Aerodynamic forces are as discussed in the section 
entitled “Finite-Difference Errors.” 

(2) To achieve a sweepback angle, the wing is rotated 
about a vortical axis through the intersection of the unswept 
elastic axis and the side of the fuselage. The tip is, however, 
terminated parallel to the airstream so that only the span 
measured along the elastic axis is unchanged in length. 

(3) Structural properties of the wing are unchanged by 
sweepback. 


Figure 15. — Flutter speed as funotion of sweepback angle 

(4) The center of mass of the fuselage is moved aft as the 
sweepback angle is increased so that it coincides roughly 
with the center of pressure of the wing. 

(5) No modifications were made for aerodynamic forces 
at the tip. 

The five basic configurations.are shown in table XII, which 
gives reference numbers for the 17 cases. The results are 
given in table XV and figure 15. Flutter characteristics of 
the two fighters show a reasonable correlation, and, in general, 
a decrease in flutter speed for sweepback angles other than 
zero. However, the bombers do not show a correlation with 
the fighters or with each other. It is significant that a sub- 
stantial change in flutter speed with sweepback angle was 
observed. In one case, flutter speed increased more than 
60 percent for a 45° sweepback, while other cases showed a 
30-percent decrease for sweepback angle of about 25°. 

In addition to the cases above it is possible to crossplot the 
variation of flutter speed with sweepback angle for the follow- 
ing parameter variations of bomber B: Wing mass density, 
wing pitching inertia, bending rigidity, torsional rigidity, 
and center-of-mass location. Most of these are plotted in 
figure 16. It is interesting to note that the general trend 
for bomber B is to a great extent independent of these varia- 
tions. 

CONCENTBATED-MASS PITCHING FLEXIBILITY 

The engines on .present-day bombers are sometimes 
mounted in nacelles on pylons some distance below the wing. 
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Figure 16. — Effect of sweepback on various modifications of bomber B 

Because of the inherent flexibility in such a structure and its 
fastening to the wing, the dynamic characteristics of the 
engine are altered. Because of the symmetry of the struc- 
ture, it is possible to write two sets of equations for the 
nacelle, one involving pitching, vertical, and fore and aft 
motion and the other involving lateral, rolling, and yawing 
motion. These sets are uncoupled except through inter- 
actions with the wing. The characteristics represented by 
the equations involving pitch have a greater effect on flutter 
characteristics, or, stated in another way, the assumption of 
a rigid pylon for lateral motion has not ordinarily been ob- 
served to introduce great differences in flutter characteristics. 
This assumption becomes less valid for wings with large 
sweepback. On the other hand, a signifi cant variation may 
be observed as the pitching flexibilities are varied. For 
pitching motion it is usually quite accurate to assume an 
effective center of rotation at some point in front of and below 
the elastic axis. Unless a specific case is to be considered, 
however, it is just as satisfactory to assume this center of 
rotation at the elastic axis, since variation in the location of 
this point has only a second-order effect. Consequently, 
in this study the center of rotation for pitching motion was 
established at the elastic axis and the pylon was assumed 
rigid for lateral motion. 

Six cases shown in table XU were investigated : 

(1) Bomber A, A=0°, mass at 0.46 span, case 176 

(2) Bomber A, A=0°, mass at tip, case 179 

(3) Bomber B, A=0°, mass at 0.46 span, ease 177 


(4) Bomber B, A=0°, mass at tip, case 180 

(5) Bomber B, A=30°, mass at 0.46 span, case 178 

(6) Bomber B, A=30°, mass at tip, case 181 

In all of these cases the chordwise position of the mass was 
basic, 60 inches forward of the elastic axis. 

In presenting the results an effort has been made to put 
the data in dimensionless form. Thus, the flutter speed is 
given as a per unit value of the flutter speed with rigid 
connection. This basic flutter speed can bo found in table 
XIV. The flexibility is conveniently measured by the 
normal-mode vibration frequency of the nacelle with the 
wing held rigid in pitch. However, instead of using the value 
of frequency in cycles per second, this frequency is measured 
in per unit value of the flutter frequency with rigid connec- 
tion. Values for the flutter frequency with rigid connection 
can also be found in table XIV. There are two frequencies 
of the nacelle which might be regarded as significant. One 
of these is the cantilever frequency in which the wing is held 
rigid in both pitch and plunge. However, for large bombers, 
the wing has such great flexibility in vertical bending that 
greater significance might be attached to the frequency when 
pitching motion is constrained and vertical motion is com- 
pletely unrestrained. Because of the location chosen for 
the basic mass, the difference in these frequencies is a factor 
of 2, the frequency with vertical motion unrestrained being 
higher. For presentation of data this higher value of fre- 
quency was chosen, because in those cases where a “tuning” 
effect was observed the maximum effect occurred when this 
frequency was equal to the flutter frequency for the basic 
rigid mass. One exception to this is observed in the discus- 
sion below. 

Results are plotted in figure 17. Nine of the twelve cases 
show a predominant decrease in flutter speed as the rigidity 
is reduced from an infinite value. Seven of these cases show 
a minimum flutter speed when the nacelle frequency is nearly 
equal to the rigid flutter frequency. This decrease varies 
between 7 and 37 percent with an average value of 18 percent. 
Two cases show a decrease in flutter speed but no tuning 
effect. The maximum rate of decrease occurs, in fact, when 
the nacelle frequency is far below the rigid flutter frequency. 
In both cases the flutter speed drops to an asymptotic value 
about six-tenths of the basic value. 

Three of the twelve cases show an increase in flutter speed 
as the rigidity is reduced from an infinite value. In two 
cases increase takes place in the region where nacelle fre- 
quency is roughly equal to the flutter frequency, and in both 
cases the flutter speed increases more than 60 percent. Tho 
last anomalous case shows a resonance or tuning effect. 
It is anomalous for two reasons: (1) The flutter speed rises 
to a sharp peak about 10 percent above basic value, and 
(2) this occurs when nacelle frequency is twice as great as 
the flutter frequency. It should be pointed out that, for 
this rigidity, the flutter frequency is equal to the nacelle 
frequency with wing attachment constrained in bending ae 
well as pitch. 

EFFECT OF A CONCENTRATED MASS 

Many aircraft structures have engines, stores, or external 
fuel tanks attached to the wing in such a way that they act 
dynamically as concentrated masses. It has long been 
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(a) Mass at 0.46 span. 

(b) Mass at tip. 

Floras 17. — Flutter speed as a function of concentrated-mass 
pitching frequency. 

known tkafc the location of such a mass has a significant 
effect on flutter. Unfortunately, other aerodynamic and 
some structural problems do not permit location of such a 
mass so that maximum flutter speed is obtained. On the 
other hand, within the restrictions imposed by other con- 
siderations, it is often possible to improve flutter character- 
istics significantly by proper choice of mass locations. 

This investigation has included a detailed examination of 
the effect of a concentrated mass on the flutter character- 
istics of several configurations of the four basic airplane 
wings. Preliminary study of this effect showed such interest- 
ing and unusual effects that the scope of the investigation 
was expanded beyond that originally proposed. The result- 
ing data are so voluminous that it is difficult to present them 
effectively. In particular, it is impractical to construct a 
table which gives all of the data obtained, and so graphical 
presentation is required. Two methods have been adopted 
in this report. For a given spanwise location of the mass, 
the flutter speed can be plotted as a function of the chordwise 
location. This has been done for all cases investigated. 
Since the concentrated mass is alined with the airstream, it 
is most convenient in cases with sweepback to move the mass 


parallel to the airstream rather than perpendicular to the 
elastic axis. Where sufficient data are available, these curves 
can be summarized in a single diagram in which lines of 
constant flutter speed are shown on a drawing of the wing 
plan form. For the concentrated mass located anywhere on 
such a contour line, the flutter speed will be the same. The 
result is essentially a topographic map of the flutter-speed 
surface, where each point on the plan form represents a 
possible location for the concentrated mass. 

Several difficulties arise with both methods of presentation. 
The main source of difficulty lies in the fact that several 
important flutter roots exist for a Aving with concentrated 
mass. For certain locations of the mass, one root will show 
lowest flutter speed, while for other locations another root 
will have the lowest flutter speed. Since the analog com- 
puter is essentially an electrical model, it is usually impossible 
to find one of these flutter speeds if another root has a flutter 
speed far below the first. It is possible, therefore, to find Avith 
certainty only those portions of a given flutter-root surface 
which he beneath all other flutter-root surfaces. For one 
configuration studied, four such distinct surfaces were posi- 
tively identified and it was not possible to establish that 
surfaces appearing at Avidely separated regions were or were 
not related. In most cases the roots were differentiated by 
obtaining essentially marginal stability for two distinct roots 
along the line where the two surfaces intersect. It can be 
readily appreciated that many points are required to establish 
the flutter-speed contours, particularly where several inter- 
secting surfaces are involved. It was, in fact, impossible in 
the time available to obtain sufficient data to establish all 
interesting features about these contours. However, it is 
believed that all important features are shown correctly in 
the figures presented here. 

The curves which shoAv flutter speed as a function of chord- 
Avise position at a fixed span frequently show intersections 
between different flutter roots. In identifying these roots, 
it is useful to know the flutter frequency associated Avith 
each root. The simplest way to present these data is to show 
the value of frequency at a few selected points along each 
curve. Where roots intersect and both frequencies were 
measured, both values are sboAvn. In some cases where 
actual- frequencies were not measured, low, medium, or high 
frequency is indicated. 

For convenience the concentrated mass was placed at the 
center of a finite-difference cell except in three cases where 
additional information was obtained by placing it halfway 
between cells. It is convenient to identify these locations 
by the cell number as has been done in table XIV, which 
assigns a case number to each configuration. It must be 
remembered, however, that the cell divisions are slightly 
different for bomber and fighter airplanes, and therefore the 
spanwise station for a given cell number Avill he different. 
The location of these stations in terms of unit span is given 
in table XLLJL. In the figures, the spanwise position is cor- 
rectly given as a fraction or per unit value of the Aving semi- 
span. 

The size of the concentrated mass and its pitching and 
rolling inertia also affect the flutter speed. Since past experi- 
ence has shoAvn that rolling inertia has a small effect, a few 
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cases were chosen for further investigation of the magnitude 
of this effect. For a concentrated mass located in the wing, 
it is reasonable to assume a radius of gyration which is a 
small fraction of the average half chord. For a mass sus- 
pended below the wing it is unlikely that the distance will 
exceed half of the average half chord. Two radii of gyration 
equal to 0.1 and 0.5 times the average chord for the fighter 
planes were chosen. In all cases considered there was no 
significant difference in flutter characteristics when the rolling 
inertia was varied from zero to the maximum value. The 
variation was, in fact, so insignificant that none of the data 
is presented in this report. In what follows it may he as- 
sumed that the rolling inertia of the concentrated mass has 
any value between the above limits. Since the mass of the 
concentrated mass and its pitching inertia have a greater 
effect, it is necessary to consider variations in these quan- 
tities in several typical cases. The basic values for mass and 
pitching inertia (or radius of gyration k) have been discussed 
in the section entitled “Characteristics of Four Represents-, 
tive Aircraft” and are given in table Xm. The variations 
of these values are summarized in table XIV. 

The most logical way to give the results is to present 
first the flutter characteristics for the basic mass on each 
particular airplane. Five airplanes were chosen: 

(1) Fighter A, A=0° 

(2) Fighter A, A=45° 

(3) Bomber A, A=0° 

(4) Bomber B, A=0° 

(5) Bomber B, A— 30° 

It was difficult to choose a typical radius of gyration for a 
mass on a fighter plane. The values used, 6, 15, and 30 
inches, are shown in table XIV. Figures 18 (a) and 18 (b) 
show the effect of chordwise location of the mass at five 
spanwise positions for fighter A with A=0°. Results for 
all three radii of gyration in pitch are plotted on the same 
sheet using different symbols for each value. Circles are 
used for the smallest value, k = 6 inches; triangles are used 
for k = 30 inches; a solid line with no symbols is used for 
the intermediate value. Abscissas for all curves are chord- 
wise distance from the elastic axis measured as per unit 
value of the wing semispan. Similar data for fighter A, 
A=45°, are shown in figures 18 (c) and 18 (d); six spanwise 
stations were used in this case. One surprising feature can 
be noted in all of these figures: The characteristics are 
relatively independent of the pitching inertia, even though 
this inertia is varied from a very large value (&=30 inches) 
to nearly zero (/c=6 inches). This does not mean that at 
any particular point the flutter speeds are identical, but the 
overall shapes of the curves are remarkably similar. 

Although these figures (figs. 18 (a) to 18 (k)) give a good 
picture of the flutter characteristics, it is easier to interpret 
the results if all data are combined to construct flutter con- 
tours as discussed earlier. Such contours for the minimum 
and maximum values of k are shown in figures 19 (a) to 
19 (h). These figures support the following conclusions: 

(1) A chordwise position aft of the elastic axis is almost 
always undesirable. 

(2) The 30- to 50-percent-span position and the tip location 
are generally undesirable. 


(3) A position forward of the elastic axis and near the 70- 
to 80-percent-span position will, in general, greatly increase 
the flutter speed. I 

Since time did not permit a complete study of the charac- 
teristics for fighter B, data were obtained only for the cases 
of 0° and 45° sweepback with mass at the tip. Comparison 
of the results shown in figure 18 (e) with the corresponding 
data for fighter A in figures 18 (a), 18 (b), 18 (c), and 18 (d) 
shows that for this location there does not seem to bo any 
significant difference in characteristics. Whether it iB safe 
to extrapolate this result to other mass locations cannot 
be said at this time. H " 



(a) Fighter A; A=0°; symmetric. 

Figure 18. — Flutter characteristics against chordwise position of con- 
centrated mass. HF, high frequency; MF, medium frequency; 
LF, low frequency. 



IN COMPRESSIBLE FLUTTER CHARACTERISTICS OF REPRESENTATIVE AIRCRAFT WINGS 


1413 


Since the study of fighter A showed that the pitch radius 
of gyration had a small effect and since the pitch radius of 
gyration of a bomber engine is relatively well defined, it 
was decided to use only one value in the study of bomber 
airplanes. However, the practice of using one engine or 
two engines on a single pylon, as well as the different sizes 
of engines, gives a possible variation in mass which might 
well exceed a factor of 2. All bomber data were therefore 
obtained with both basic mass and half basic mass. For 
ease of comparison, the two sets of -data are plotted side by 
side in the figures. Data for bomber A, A=0°, are contained 
in figures 18 (f) and 18 (g). Results for bomber B, A=0°, 
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(b) Fighter A; A=0°; antisymmetric. 
Fiqtjbb 18. — Continued. 
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(d) Fighter A; A=45°; antisymmetric. 
Figure 18. — Continued. 


figures shows some deviations from the results for fighter 
airplanes. The following conclusions can be drawn: 

(1) A position aft of the elastic axis is almost always 
undesirable. 

(2) The tip region is generally undesirable as a location 
for the mass. 

(3) With few exceptions, any position forward of the 
elastic axis and between the root and the 90-percent-span 
position will give flutter speed equal to or greater than the 
bare-wing flutter speed. 


(4) There are, in most cases, no practical locations which 
give any significant improvement in flutter characteristics. 
Two cases will be noted in which the speed might be increased 
40 percent. The others are restricted to a 10- or 20-percent 
improvement. 

Since fighter planes showed remarkable variation in 
flutter characteristics with mass position, it was believed 
to be necessary to examine the effect of changes in the size 
(mass) of the concentrated mass. This was first studied at 
two spanwise positions for fighter A, A=0°. The positions 
are the tip and station 5 (0.79 span) . Flutter characteristics 
as functions of chordwise position were measured for several 
values of mass. The three values of radius of gyration 
given in table XD1 were used for all mass values, oxcept 
that for very small values of mass only the 6-inch value was 




(d) Continued. 

(e) Fighter B; data taken at tip. 
Figurb 18. — Continued. 
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(f) Bomber A; A=0°; symmetric. 
Figtxbe 18. — Continued. 



(g) Bomber A; A=0°; antisymmetric. 
Figure 18. — Continued. 
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(h) Bomber B; A=0°; symmetric. 
Figure 18. — Continued. 
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Chordwise position 


(j) Bomber B; A=30°; symmetric. 
Figube 18. — Continued. 


(k) Bomber B; A=30°; antisymmetric. 
Figttkb 18. — Concluded. 
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(a) Fighter A; A=0°; symmetric; radius of gyration, 6 inches; 

1.0t>/=989 miles per hour. 

(b) Fighter A; A=0°; symmetric; radius of gyration, 30 inches; 

1.0»/=989 miles per hour. 

(c) Fighter A; A=0°; antisymmetric; radius of gyration, 6 inches; 
1.0»/= 1,091 miles per hour. 

* Figure 19. — Contours of constant flutter speed with concentrated 

mass. 

used. However, because of the similarity of results, data 
are 'presented only for the maximum value (£=30 inches) 
andjminimum value (£= 6 inches). Case numbers are listed 
in table XIV. 

The results for tip location shown in figure 20 (a) show a 
very interesting progression in characteristics as the mass is 



(d) Fighter A; A=0°; antisymmetrio; radius of gyration, 30 inches; 

1.0t>/= 1,091 miles per hour. 

(e) Fighter A; A=45°; symmetric; radius of gyration, 6 inches; 

1.0t>/= 1,028 miles per hour. 

Figure 19. — Continued. 

reduced to zero. Most striking is the fact that no significant 
change takes place when the mass is varied from twice basic 
value to half basic value. Even with mass reduced to 8 per- 
cent of basic value, the three flutter roots for antisymmotric 
motion and the two flutter roots for symmetric motion 
can still be identified though their characteristics are by this 
time somewhat altered. Similar data for the mass at 0.79 
span are presented in figure 20 (b). Figure 20 (a) also 
illustrates graphically the danger in extrapolating results. 
For a mass at the 0.10 chordwise position and symmetric 
motion, a mass equal to 0.008 basic mass increases flutter 
speed to 1.07. Doubling the mass increases it to 1.16. 
Again doubling it will increase the speed to some unknown 
value greater than 1.40. However, if the mass is again 














(f) Fighter A; A =45°; symmetric; radius of gyration, 30 inches 

1.0»/= 1,028 miles per hour. 

(g) Fighter A; A=45°; antisymmetric; radius of gyration, 6 inches 

1.0»/=881 miles per hour. 

FraxjBB 19. — Continued. 
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(k) Bomber B; symmetric; 1.0t>/= 1,216 miles per hour. 

(l) Bomber B; antisymmetric; 1.0 p/= 1,216 miles per hour. 

(m) Bomber B; symmetric; 1.0t»/= 1,165 miles per hour. 

(n) Bomber B; antisymmetric;. 1.0t>/= 1,244 miles per hour. 

Figure 19. — Concluded. 



(a) Concentrated mass at tip. 

Figure 20. — Effect of size of concentrated mass on flutter character- 
istics. Fighter A; A=0°. 


increased by a factor of about 2.8, the flutter speed (of 
another root) will have dropped to 1.00 again. 

Flutter characteristics were also measured for basic and 
half basic mass at the wing tip with sweepback angle of 46°. 
This case was chosen because the antisymmetric character- 
istic for a 30-inch radius of gyration (fig. 18 (d)) showed a 


INCOMPRESSIBLE FLUTTER CHARACTERISTICS OF REPRESENTATIVE AIRCRAFT WINGS 


1421 



(b) Concentrated mass at 0.79-span position. 
Figure 20. — Concluded. 


very unusual characteristic. Figure 21 (a) shows that reduc- 
tion of the mass by a factor of 2 eliminates the anomalous 
behavior, but in all other respects gives results which are 
essentially the same as the basic mass. Figure 21 (b) gives 
similar data for the mass at the 0.79 span position with A= 
45°. Again the results for basic mass and half basic mass are 
not significantly different. It is perhaps unwise to attempt 
any general statement, but there is every indication that the 
essential features shown in the flutter-speed contours of 
figure 19 would not be profoundly altered if either mass or 
pitching inertia were increased or decreased by a factor as 
large as 2.0. 

ACCURACY AND PROBABLE ERRORS 

It is behoved that a discussion of the expected accuracy 
of the results will be more meaningful now that the reader has 
observed the nature and scope of the data obtained. The 
analog computer is not composed of perfect electrical 
elements. For example, the inductors used in this study 
have loss characteristics corresponding to a damping factor 
of about {7=0.01. Transformers also have significant losses.. 
The electrical analog of the airplanes studied in this report 
had an electrical damping corresponding to a structural 
damping between {7=0.02 and {7=0.03. This is not greatly 
different from the damping to be found in conventional 
aircraft construction, so no corrections were made for this 
internal damping. 

No general statement can be made about the effect of 
random computer errors. Some give rise primarily to an 
error in the damping factor of the roots, in which the [g—v) 
curve is shifted vertically. Other errors give rise basically 
to an error in velocity, in which the curve is primarily 
shifted horizontally. Since the slope of the curve of g 
against v is by no means constant, it becomes impossible 
to give a specific figure for accuracy of flutter speed. In 
unusual cases, where roots are of the type shown in figure 12, 
flutter may be predicted when in fact it will not occur for 
these roots at all. This distinction, which must be made 
mathematically, is of no importance in practical cases. An 
airplane which shows a flutter damping of 0=0.01 and is 
therefore theoretically stable is not to be regarded as any- 
more satisfactory or useful than one which shows a damping 


0=— 0.01 and would therefore theoretically fly apart. 
Manufacturing tolerances and the safety factors required in 
aircraft will not permit use of an aircraft unless it is moder- 
ately stable for a significant variation in all structural 
parameters. It is believed that, exclusive of errors intro- 
duced by the finite-difference structure and approximations 
in the aerodynamic theory, the results obtained in this study 
have a probable error in damping factor of about 0= ±0.02 
or a probable error in flutter speed of about 2 percent, which- 
ever is applicable in the light of the above discussion. How- 
ever, trends obtained by variation of structural parameters 
are considerably more accurate than this would imply, 
since any error would persist with roughly the same value in 
all cases involving such parameter variations. 

For cases 1 to 181, it is possible to construct curves of g 
against v although they have not been made a part of this 
report. These curves serve to show the steepness with which 
the roots pass through flutter and give some indication of the 
accuracy of the flutter speed. This situation does not exist 
for cases 181 to 289. For concentrated-mass variations, 



(a) Basic mass and 0.5 basic mass at tip. 

(b) Basic mass and 0.5 basic mass at 0.79-span position. 

Figure 21. — Comparison of characteristics for basic mass and 0.5 
basic mass. Fighter A; A=45°. 
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computations were carried out in such a way that only flutter 
speeds and frequency were obtained. Consequently, it is 
not possible to delineate areas which. are “safe” from the 
standpoint of flutter. It is known, for example, that where 
a long pendant lobe is observed, as for case 235 in figure 
18 (d), the system is barely unstable everywhere within this 
Lobe. There are s imil ar regions where the system is barely 
unstable within an elliptical-shaped boundary, for example, 
case 212 (fig. 20 (b)). Similarly, case 214 is barely stable 
in this region and yet no flutter root is even shown since 
the system does not become actually unstable at any point. 

These remarks are not made to show the flutter curves 
to be valueless, but to caution the reader against making 
inferences not contained in the report and not legitimately 
supported by the data presented here. 

SUMMARY OF RESULTS 

Certain trends in flutter characteristics for typical modem 
aircraft seem to be indicated by-this study of the incompres- 
sible flutter characteristics of aircraft wings. In some cases 
a few deviations are found. It is probable that if more 
extreme aircraft designs were considered even more would be 
observed. Nevertheless, this summary may serve as a 
useful guide. 

(1) In the following table are fisted the average changes in 
flutter speed for a 1-percent change in each parameter, the 
change being made from the basic value except for center- 
of-mass and elastic-axis location, in which cases several 
locations are assumed. 


Change in parameter 

Average change 
in flutter speed, 
percent 

Increase wing mass 1 percent - - 

0 


— . 1 

Increase bending rigidity 1 percent 

0 


. 5 

Increase wing chord. 1 percent 

— . 7 

Center of mass forward 1 percent of chord 
from — 

40-percent location _ _ _ 

3. 0 

50-percent location- 

1. 7 

60-percent location 

. 8 

Elastic axis forward 1 percent of chord from — 
40-percent location - 

3. 2 


1 . 8 



(2) Localized change in torsional rigidity is most effective 
in changing flutter speed of a bare-wing airplane if the 
change is made between midspan and root. For bombers 
with large concentrated mass on the wing, torsional rigidity 
either inboard or outboard of the mass will usually govern 
flutter speed, depending on the type of flutter existing. 
Sweepback was not observed to have a systematic effect. 
Pitching flexibility of the concentrated-mass support has a 
definite influence on flutter speed. In many cases a tuning 
effect was observed, with a 10- to 40-percent decrease in 
flutter speed. This effect was not always observed; in 
some cases, the flutter speed was significantly increased. 


(3) Perhaps the most interesting results will be found in 
the effects of a concentrated-mass location. For a wide 
range of mass and radius-of-gyration values the results 
were very systematic. For fighter-type planes it was found 
that: 

(a) Aft chordwise positions are usually undesirable. 

(b) The 30- to 50-percent-span and tip locations are 

generally undesirable. 

(c) A forward location near the 70- to 80-percent-span 

position will, in general, greatly increase flutter 
speed. 

For bomber-type planes these results are somewhat modified: 

(a) Aft chordwise positions are usually undesirable. 

(b) The tip location is generally undesirable. 

(c) "With few exceptions any location forward of the 

elastic axis and between root and 0.90-span posi- 
tion is satisfactory, although flutter speed is 
rarely greater than bare-wing flutter speed by 
any significant amount. 

There is perhaps no need to remark that these results can 
be altered by introduction of a flexibility in the concentrated- 
mass support. 

California Institute of Technology, 

Pasadena, Calif., June 6, 1965. 
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